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ABSTRACT 


The  net  influence  of  sab-grid  or  smaller  unresolv&ble  scale 
processes  (which  cannot  be  directly  detected  from  conventional 
observations)  on  the  tangential  momentum  field  in  tropical  cyclones  has 
been  calculated  as  a  residual  from  the  grid-scale  momentum  budgets  of 
sixteen  rawinsonde  composite  data  sets.  Results  show  that  cyclone 
systems  of  the  Pacific  and  Atlantic  that  are  undergoing  development  have 
similar  residual  profiles  in  the  vertical:  positive  cyclonic 
acceleration  at  upper  and  lower  levels  and  negative  cyclonic 
acceleration  at  between  levels.  Stronger  systems  have  a  greater 
magnitude  of  the  sub-grid  scale  effect.  These  results  are  believed  to 
represent  primarily  the  cumulus  convective  transport  (defined  as  cumulus 
friction).  A  simple  cloud  model  is  developed  to  diagnose  momentum 
rearrangement  in  the  vertical  through  mass  recycling.  Model  results 
agree  well  with  the  residual  calculations  from  the  composites. 

The  total  convective  effect  (TCE)  on  the  tangential  momentum 
includes  the  vertical  transports  by  all  the  upward  and  downward  motions 
which  can  be  decomposed  into  the  mean  vertical  motion  and  the  local  mass 
recycling.  This  TCE  is  found  to  reduce  both  the  low  level  cyclonic  and 
the  upper  level  anticyclonic  winds.  These  results  imply  that  the  net 
effect  of  the  cumulus  convective  transport  is  down-gradient  and  acts  to 
reduce  the  vertical  wind  shear  of  the  cyclone.  This  mutual  interaction 
between  the  cumulus  convection  and  the  larger  scale  flow  features  and 
their  possible  influence  on  tropical  cyclone  development  is  discussed. 
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1 .  INTRODUCTION 


The  importance  of  cumulus  convective  processes  on  the  horizontal 
momentum  fields  of  tropical  weather  systems  has  been  discussed  by  Gray 
(1967,  1973),  Reed  and  Johnson  (1974),  Steven  £jt  a_l.  (1977),  Cho  et_a_l*. 
(1979)  and  Shapiro  and  Stevens  (1980).  It  is  impossible  to  directly 
determine  such  convective  influences  because  cumulus  convection  occurs 


on  spa 


ce  and  time  scales  significantly  smaller  than  those  in  which 


measurements  of  the  large-scale  flow  (defined  here  as  that  determined  by 
a  conventional  rawinsonde  observational  network.)  patterns  are  taken.  As 
a  result,  many  investigators  have  tried  to  determine  the  residual  of 
the  large-scale  momentum  and  vorticity  budgets  (Williams  and  Gray,  1973, 
Reed  and  Johnson,  1974;  Shapiro,  1978;  Cho  irt  aj.. ,  1979  and  Stevens, 
1979).  If  cumulus  convection  is  the  only  important  small  scale  process, 
this  residual  may  be  considered  to  be  largely  contributed  by  the  cumulus 
convective  transport. 

Theoretical  cloud  models  have  also  been  developed  to  diagnose  the 
influence  of  cumulus  convective  transports  on  the  large-scale  dynamic 
fields  (Austin  and  House,  1973;  Feed  and  Johnson,  1974;  Shapiro,  1978; 
Cho  et  al.»  1979;  Shapiro  and  Stevens,  1980).  The  model-diagnosed 
results  are  often  compared  with  the  calculated  residuals  of  the  large- 
scale  budgets.  A  severe  inconsistency  between  the  model-diagnosed 
results  and  calculated  residuals  would  throw  doubt  on  any  conclusions 
based  on  either  the  calculated  residuals  and/or  model-diagnosed  results. 
However,  a  reasonable  consistency  between  these  two  would  add  support  to 
both  observational  and  model-inferred  results  and  help  support  the 
hypothesis  that  cumulus  convection  is  the  most  important  small  scale 


process.  Such  results  would  also  demonstrate  that  the  calculated 


residual  can  well  approximate  the  cumulus  convective  transport.  These 
cloud  models  might  also  provide  a  framework  for  a  parameterization 
scheme  to  represent  cumulus  convective  transport  in  numerical  models 
(e.g.  Stevens  ejt  a_l . ,  1977;  Fingerhut,  1982;  Hack  and  Schubert,  1980). 

Cumulus  convection  is  typically  observed  to  be  quite  substantial  at 
the  trough  region  of  a  tropical  synoptic  wave  and  at  inner  radii  of 
tropical  cyclone  systems.  Cumulus  convective  transports  should 
therefore  be  both  detectable  in  and  important  to  the  dynamics  of  these 
systems.  Unfortunately,  these  weather  systems  typically  exist  over 
tropical  oceans  where  conventional  observations  are  limited.  It  is 
usually  not  possible  to  determine  quantitatively  the  large-scale  flow 
patterns  (particularly  the  large-scale  mean  vertical  motion)  in 
individual  case  situations.  One  is  obliged  to  composite  many  similar 
cases  together  to  obtain  quantitative  and  meaningful  large-scale 
circulation  patterns. 

This  composite  technique  has  been  used  by  Reed  and  Recker  (1971)  to 
study  tropical  synoptic  waves  over  the  western  Pacific  using  3  months  of 
rawinsonde  data.  The  data  of  the  GARP  (Global  Atmospheric  Research 
Project)  Atlantic  Tropical  Experiment  (GATE)  have  also  been  composited 
in  the  studies  of  tropical  synoptic  wave  systems  over  Africa  (Reed  ej: 
al .  .  1977)  and  the  eastern  Atlantic  Ocean  (Thompson  e_t  a_l. ,  1979).  The 
rawinsonde  compositing  at  Colorado  State  University  (CSU)  has  produced 
quite  meaningful  quantitative  analyses  of  large-scale  horizontal  and 
vertical  motion  features  of  tropical  cloud  clusters  and  tropical  cyclone 
systems  (Williams  and  Gray,  1973;  Frank,  1977a,  1977b;  McBride,  1981; 
McBride  and  Zehr,  1981;  and  other  papers). 
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Assuming  that  such  composited  large-scale  motions  are  realistic, 
one  may  proceed  a  step  further  and  ask  the  question:  to  what  extent  «f* 
the  small  (or  unre sol vable )  scales  of  motion  an  important  component  to 
the  necessary  quantitative  balances  of  the  large  (or  resolvable)  scales? 
To  what  extent  can  the  large-scale  momentum  (also  energy  and  moisture) 
budget  be  made  without  considering  the  smaller  scale  convective 
processes?  The  answers  to  these  questions  are  of  fundamental  importance 
to  the  understanding  and  modeling  of  such  convective  weather  systems. 

Many  efforts  have  been  made  to  answer  such  questions  for  tropical 
waves  using  composite  data  (Reed  and  Johns  on.  1974;  Shapiro.  1978; 
Stevens,  1979;  Cho  e_t  a_l .  ,  1979  and  Shapiro  and  Stevens,  1980). 

However,  no  one  has  yet  performed  such  a  study  for  the  region  of  the 
tropical  cyclone  beyond  the  eye-wall.  This  paper  thus  attempts  to 
answer  the  above  questions  for  the  inner  1-5°  latitude  radius  region  of 
tropical  cyclone  systems  -  including  tropical  cloud  clusters  and 
hurricanes  and  typhoons. 

To  determine  the  cumulus  momentum  transport  as  a  residual  of  the 
large-scale  momentum  budget,  one  is  forced  to  evaluate  the  pressure 
gradient  force  which  is  usually  less  reliable.  Therefore,  most  previous 
studies  (mentioned  above)  have  resorted  to  an  analysis  of  the  vorticity 
budget  in  which  the  pressure  gradient  force  is  not  involved.  However, 
these  vorticity  budget  calculations  include  a  twisting  term  which 
requires  a  knowledge  of  the  horizontal  gradient  of  vertical  motion. 

This  twisting  term  was  neglected  in  the  analysis  of  Reed  and  Johnson 
(1974).  Shapiro  (1978)  calculated  this  term  and  suggested  that  this 
twisting  effect  is  perhaps  not  important.  Stevens  (1979)  also  found 
that  the  twisting  term  was  smaller  than  the  other  terms  in  the  large- 
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seals  vorticity  budget  equation.  On  the  other  hand,  Cho  e£  al .  (1979) 
and  aore  recently  Yanai  (private  coaaunication)  have  found  that  the 
twisting  tern  is  often  very  iaportant  in  their  analyses  of  tropical 
conveotion.  These  apparently  contradicting  conclusions  indioate  soae 
general  doubt  in  the  reliability  of  the  vorticity  analysis.  Such 
coaplication,  however,  will  not  arise  in  the  present  study  which  uses  a 
cylindrical  rather  than  a  rectangular  coordinate  systea. 

In  this  study,  the  tangential  aoaentua  budget  for  the  inner  region 
o 

(1-5  latitude  radius)  of  tropical  cyclone  systeas  is  considered  rather 
than  the  vorticity  equation.  After  taklna  an  ariauthal  averaae  of  the 
larae-scale  tanaential  wind  equation,  no  knowledte  of  the  heiaht  or 
pressure  field  is  required  because  the  pressure  aradient  tera  drops  out 
in  an  intearation  around  the  vortex.  No  complicated  twisting  term  is 
involved  in  the  calculation  of  the  residual. 

The  calculated  residual,  or  the  unresolvable  part  of  the  total 
vertical  transport  process,  is  defined  as  cuaulus  friction,  denoted  by 
cuFg  (after  Schneider  and  Lindzen,  1976;  Stevens  et  al.  1977).  Values 
of  this  cuaulus  friction  for  16  coaposited  tropical  cyclone  data  sets 
are  presented  in  Chapter  3.  The  consistencies  and  systeaatic  variations 
between  data  sets  are  also  discussed.  In  Chapter  4  these  results  will 
be  coapared  with  those  diagnosed  froa  s  staple,  single  cloud  aodel. 

The  total  convective  effect  (TCE),  which  includes  both  the 
resolvable  and  the  unresolvable  parts  of  the  verticsl  transport  of 
tangential  aoaentua,  is  also  determined.  Results  are  presented  in 
Chapter  5.  The  autual  influences  of  the  sasll  scale  convective  systeas 
and  the  large-scale  aoaentua  fields  sre  also  discussed. 


5 


2.  RAW IN SONDE  COMPOSITE  DATA  SET 


The  rawinsonde  composite  data  sets  used  in  this  study  were 
established  by  Prof.  V.  M.  Gray's  tropical  cyclone  research  group.  They 
consist  of  14  years  (1961-1974)  of  northwest  Atlantic  and  10  years 
(1961-1970)  of  northwest  Pacific  rawinsonde  soundings.  Soundings  are 
composited  in  a  cylindrical  coordinate  system  whose  center  is  fixed  on 
the  cyclone  center.  Data  are  resolved  into  7  radial.  8  azimuthal  and  17 
vertical  grid  points.  Each  grid  point  value  on  each  level  represents  an 
area  average  over  2°  latitude  radius  and  45°  azimuth  as  shown  in  Fig.  1. 


Fig.  1.  Cylindrical  grid  for  composite  data.0  Arrow  points  toward 
north.  The  radial  grid  spacing  is  2  latitude.  Bach  grid 
point  value  represents  the  average  of  all  soundings  in  the  grid 
box  (shaded  area). 


The  16  composite  tropical  cyolones  used  in  this  study  are  summarised  in 
Table  1.  For  a  more  complete  description  of  these  data  sets,  the 
Colorado  State  University  (CSU)  rawinsonde  composite  philosophy,  and 


Table  1  (continued) 
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4.  mature  typhoon  or  hurricane  5,  super typhoon 


8 


specific  composite  techniques,  the  render  is  referred  to  Villisas  end 
Gray  (1973),  Frank  (1977s),  MoBride  (1981),  Gray  (1981)  and  Gray  et  si. 
(1982). 

2.1  Smoothini  of  the  Composite  Data 

The  composite  data  fields  are  considered  to  represent  the  average 
state  of  a  tropical  cyclone.  However,  it  is  inevitable  that  the 
composite  average  do  include  a  degree  of  meteorological  noise.  To 
reduce  this  noise  the  smoothing  rontines  contained  in  the  SPLPAJf  package 
available  from  the  National  Center  for  Atmospheric  Research  (NCAR)  vere 
need.  This  package  contains  rontines  for  fitting  (by  least-sqnares)  a 
multi-dimensional  cubic  spline  to  values  at  specific  points,  and  those 
for  evaluating  the  spline  at  any  point.  A  2-dimensional  (z  and  y- 
directions)  cubic  spline  was  used  to  smooth  the  composite  data  at  each 
level.  The  degree  of  'smoothness'  depends  on  the  number  of  nodes  in 
each  direction  -  the  more  nodes,  the  less  smoothed  the  data.  In  order 
to  choose  an  optimum  degree  of  smoothness,  the  azimuthal  average  of  the 
original  composite  tangential  wind  at  900  mb  and  2S0  mb  for  the  WPD4 
data  set  and  the  smoothed  tangential  winds  using  5-node  and  9-node 
smoothing  schemes  were  compared.  The  results  are  shown  in  Fig.  2.  The 
five-node  smoothing  scheme  used  by  Nunez  (1982)  is  sdequate  for  outer 
radii  (defined  as  radius  >  5°  latitude)  where  the  radial  gradient  is 
relatively  weak;  however,  it  is  inadequate  to  resolve  radial  gradients 
at  inner  radii  (defined  as  1-5°  latitude  radius).  The  9-node  smoothed 
data  (used  in  this  study)  resemble  the  originsl  data  much  more  than  the 
5-node  smoothed  data  at  inner  radii,  which  is  the  ares  of  interest. 

One  of  the  advantages  of  using  smoothing  routines  is  the 
computstion  of  horizontal  gradient.  Another  advantage  is  that  the  value 
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a v  any  point  inside  the  smoothing  domain  can  be  interpolated.  This 
provides  a  way  of  calculating  the  tangential  wind  change  at  each  grid 
point  due  to  storm  motion.  This  is  discussed  in  more  detail  in  section 
3.1.  In  addition,  values  at  grid  points  with  missing  data  can  be 
objectively  interpolated  with  a  minimum  of  smoothing. 


Fig.  2.  Original  (solid),  9-node  (dashed)  and  5-node  (dotted)  amoothed 
tangential  wind  field  at  900  mb  and  250  mb  for  WPD4  data  set. 
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3.  CUMULUS  TRAN SPORTS  OF  TANGENTIAL  MOMENTUM  (CUMULUS  FRICTION)  CALCULATED 
AS  A  RESIDUAL  OF  THE  GRID  SCALE  BUDGET 

Without  the  exist  •no*  of  omuls*  conveotive  systems,  the  mail 
•oala  processes  of  the  tropical  cyclone  would  likely  be  of  little 
importance  compared  to  the  large  scale  processes.  However,  in  a 
tropical  cyclone  system  with  a  horizontal  scale  of  several  hundred 
kilometers,  there  are  many  convective  elements  of  horizontal  scale  of 
several  kilometers  or  larger,  which  Arnold  (1977)  referred  to  as  Basic 
Convective  Elements  (BCE).  Figure  3  shows  the  BCE's  associated  with  a 
mature  tropical  cyclone  system.  These  convective  elements  often  exist 
as  line  convection  and  frequently  propagate  at  a  speed  greater  than  the 
environmental  wind  speed  (Fujita  and  Black,  1970;  Zipser,  1977).  These 
convective  systems  transport  heat,  momentum  and  moisture  vertically  and 
horizontally  in  a  manner  which  can  be  quite  different  from  that  of  the 
larger,  resolvable  scale  circulations.  Therefore,  large-scale  dynamic 
and  thermodynamic  budgets  computed  only  from  the  transports  on  the  large 
or  resolvable  scale  will  probably  be  inadequate. 

This  paper  attempts  to  determine  how  these  small  soale  convective 
systems  can  affect  the  large  scale  momentum  field  in  a  tropical  cyclone 
system.  Two  scales  are  thus  defined:  the  grid  (or  large)  scale  and  the 
sub-grid  (or  small)  scale.  Grid  scale  processes  are  those  of  horizontal 
scale  equal  to  or  greater  than  twice  the  grid  spacing.  Sub-grid  scale 
processes  have  a  horizontal  scale  smaller  than  twice  the  grid  spacing. 

These  sub-grid  scale  influences  cannot  be  explicitly  defined  by  our 
composite  rawinsonde  data  but  they  can  be  calculated  as  a  residual  from 
the  grid  scale  budgets. 


oncnt  : 
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dt  rdr  rd0  dp  rd0 
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where  brackets  [  ]  denote  the  grid  point  quantity  which  is  implicitly 
the  grid  box  average.  The  tern  [F^]  is  defined  as  the  tangential 
friction  due  to  all  the  sub-grid  scale  effects,  or 
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_  JjgJxU  _  artu'vM  _  aiTVl  _  auvi  .  .  .  [3] 

r  rdr  rdO  dp  p 


where  priae  (')  denotes  the  deviation  froa  the  grid  box  average  and  x 

P 

is  the  sub-grid  scale  pressure  gradient  force,  which  is  probably  a 

result  of  convection-induced  dynaaic  pressure. 

The  total  vertical  transport  of  tangential  aoaentua  associated  with 

all  the  convective  processes  includes  transport  by  the  mean  vertical 

notion,  (-  ) ,  which  can  be  deterained  directly  froa  the 

op 

dfw'v'l 

composite  data  and  that  by  vertical  eddy  processes,  -  ^  ,  which 

cannot  be  directly  deterained  froa  the  conposite  data.  This  eddy  tern, 
d  Tu'v* 1 

-  .  ,  should  be  iaportant  in  the  presence  of  convective  systeas 

op 

and  is  defined  as  cuaulus  friction,  denoted  as  [CuFg],  following  the 

noaenclature  of  Schneider  and  Lindzen,  1976;  Stevens  Jl. ,  1977;  and 

Fingerhut,  1982.  However,  this  eddy  tern  is  not  exclusively  the  result 

of  cuaulus  convection.  Any  horizontal  variation  of  vertical  notion  in  a 

grid  box  is  also  included  in  this  eddy  tern. 

The  most  iaportant  sub-grid  scale  effect  is  often  thought  to  be 

associated  with  the  vertical  eddy  prooess,  -[  ,  ] .  How  the  sub-grid 

op 

scale  pressure  gradient  force,  x  ,  can  affect  the  large  scale  aoaentua 

P 

fields  is  still  not  well  understood.  However,  this  pressure  gradient 


force  should  act  in  all  directions  and  average  ont  over  a  large  area, 
and  is  thus  neglected  in  this  study  (as  in  other  studies) .  The  first 
three  terns  on  the  right  hand  side  of  Eq.  3  are  associated  with 
horizontal  snb-grid  scale  eddies,  snch  as  transient  eddies,  which  nay  be 
of  significant  nagnitude.  This  transient  eddy  effect  will  be  discussed 
later. 

Rearranging  Eq.  2  gives: 
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Except  for  the  local  tine  derivative  of  tangential  wind  and  tangential 
surface  friction,  (  and  iSFgl)  *  *H  the  other  terns  on  the  right 

hand  side  of  Eq.  4  can  be  evaluated  by  applying  finite  differencing  on 
the  conposite  data. 

3.1.1  Local  Change  of  Tangential  Monentun 

d[v] 

The  looal  tine  derivative  ^ ’  in  a  geographically-fixed 
coordinate  systen  consists  of  two  parts:  one  due  to  the  notion  of  the 
systen,  denoted  as  (3v/dt>B,  and  the  other  due  to  the  intensity  change 
of  the  systen,  denoted  as  (dv/dt)j.  Assuming  that  the  systen  is  in  a 
steady  state  and  noving  with  a  naan  speed  (c)  in  a  constant  direction, 
the  notion  part  can  be  rewritten  as: 


( 


at' 


c 


lx  .  !s! 

da  r 


where  d/8s  is  the  gradient  along  the  direction  of  notion  of  the  systen 
and  Cg  the  tangential  component  of  the  cyclone  speed.  The  right  hand 
side  of  this  equation  is  the  tangential  conponent  of  (c  *V  )V  ,  where 


the  arrow  — )  represents  a  vector  quantity.  Snoothing  routines  have 
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bean  used  to  Interpolate  the  value  of  v  at  pointa  inaide  the  composite 
grid  where  the  data  are  not  defined.  Thia  Hake a  it  poaaible  to 
oalonlate  dv/ds  at  any  grid  point  by  naing  a  finite  differenoe  scheme: 

is.  .  v(A»)  -  v(-A») 
ds  2Aa 

where  v(As)  ia  the  tangential  wind  speed  at  the  point  As  distance  ahead 
of  the  grid  point,  along  the  wowing  direction,  and  v(-As)  ia  the  point 
behind. 

To  oalcnlate  the  part  dne  to  intensity  change  of  the  syatem, 
(dv/dt)j,  the  syatem  is  assnmed  to  change  its  intensity  to  a  later 
stage  (see  Table  2)  in  a  certain  time  period.  At.  These  time  periods 
are  determined  by  referring  to  those  used  by  McBride  (1979)  and  the 
typical  rates  of  change  of  the  maximum  wind  speed  in  tropical  cyclones 
(Neal  and  Holland  jjf.  al. .  1978).  Table  2  lists  the  time  periods  (At) 
between  different  intensity  stages  of  a  cyclone.  The  last  line  lists 
those  storm  systems  assumed  to  be  in  steady  state. 

TABLE  2 

Time  period  (At)  for  a  system  to  intensify  or  weaken  from  an  original 
stage  to  the  later  stage.  Last  line  lists  those  data  sets  considered  to 
have  no  intensity  change. 


Orlainal  State 

At  (davs) 

Lilt*  S1484 

ATD1 

1.5 

ATD2 

WPD2;  ATD2 

2.0 

WPD3 ; ATD3 

WPD3;  ATD3 

3.0 

WPD4;  ATD4 

WPD4(D) ,  WPD4(F);  ATD4(D) 

1.5 

WPD4;  ATD4 

*PD3(D);  AIN 3 

2.0 

WPD4;ATN1 

WPN1,  WPN3 ,  WPD4,  WPD5,  No  Intensity 

ATN1,  ATD4  Change 
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3.1.2  Surface  Friction 

Surface  friction  needs  to  be  considered  carefully  in  a  momentum 
budget  study  and  requires  a  knowledge  of  the  surface  stress  and  its 
vertical  distribution.  Many  empirical  formulae  have  been  developed  to 
calculate  tbe  surface  stress  (Garratt,  1977).  However,  none  of  theae 
are  applicable  to  the  turbulent  boundary  layer  of  a  tropical  cyclone 
because  of  the  unknown  convective  effect  on  the  depth  of  the  aixing 
layer  and  the  breaking  ocean  wave  condition  associated  with  strong  and 
gusty  winds.  However,  over  the  1-5°  latitude  radius  domain  considered 
here  these  effects  are  assumed  to  be  of  secondary  importance.  Thua, 
following  Frank  (1977b)  the  Deacon  bulk  aerodynamic  formula  (Boll,  1965) 
is  used  to  determine  surface  stress: 

ro  "  ""Si  VloKo 

where  CD  =  1.1  x  10  3  +  4.0  x  10  5  lVi()l  ^iq  in  lv10l  T10 

are  the  total  wind  speed  and  tangential  wind  speed  at  10  m  height 
respectively.  The  wind  speed  at  900  mb  is  scaled  by  0.8  to  approximate 
the  wind  speed  at  10  m  height.  In  addition,  a  constant  multiplier,  1.1 
for  Pacific  data  sets  and  1.3  for  Atlantic  data  sets,  is  used  to  adjust 
the  surface  stress  calculated  by  the  bulk  formula,  to  account  for  the 
wind  variation  due  to  terrain  effects  which  vary  between  the  ocean 
basins.  (There  are  more  mountains  around  the  Atlantic  Ocean  basin.) 

The  surfaoe  stress  is  assumed  to  vanish  at  900  mb,  which  only 
allows  the  surfaoe  friotion  to  be  below  this  layer.  Also,  the  surface 
friction  is  assumed  to  be  larger  at  the  surfaoe  and  990  mb  than  at  900 
mb.  However,  no  matter  how  the  effect  of  the  surfaoe  friction  is 
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vertically  distributed,  its  vertical  integration  (surface  stress)  is  the 
saae,  that  is, 

l  SFe  4*  ■  '  S^oKo 

where  h  is  the  depth  of  the  nixing  layer  and  -  0.  (Density  p  is 
fairly  constant  in  the  nixing  layer.) 

3.2  Axinnthal  Averaging  of  Grid-scale  Tangential  Monentnn  Budget 

The  sub-grid  scale  eddy  effect  [Fg]  at  each  grid  point  (shown  in 
Fig.  1)  can  be  calculated  fron  Eq.  4  using  conposite  rawinsonde  data. 
However,  these  calculated  residuals  night  have  a  nagnitude  of  up  to  SO  n 
s  *day  1  especially  at  upper  levels.  These  large  and  unrealistic  values 
are  prinarily  a  result  of  the  pressure  gradient  tern  in  Eq.  4.  The 
nonentun  budget  study  by  Shapiro  and  Stevens  (1980)  had  the  sane 
problen.  This  is  the  prinary  reason  why  other  researchers  have  resorted 
to  the  use  of  the  vorticity  rather  than  the  nonentun  equation.  In  & 
cylindrical  coordinate  system,  however,  the  pressure  aradient  force  tern 
can  be  ellninated  by  takinx  the  azimuthal  averaae.  When  such  azimuthal 
averaging  is  performed,  Eq.  4  then  becomes 

(fTj  «  iM  +  Lrl  +  9£*Uy1  +  +  lallzl  .  jgpj  151 

V  vt  TOT  O  p  T  0 

where  the  overbar  (  )  denotes  an  azimuthal  average.  Note  that  this 

equation  requires  only  observed  wind  data  which  are  considered  to  be 
much  more  reliable  than  observed  height  data. 

Taking  azimuthal  averages  of  Eq.  3  gives 

to  .  _  sis .  insm  .  nzm 

0  r  rdr  dp 
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or 


2fc,v>]  at  uV]  dlm'y'l 
L0J  *  r  3r  “  3p 


[6] 


where  the  sub-grid  scale  pressure  gradient  ([x  ])  is  also  averaged  out. 

P 


Transient  eddy  effects,  defined  as 


_  [nlylJ 


,  are  estimated  by 


r 


S-X)  _  ifiltil) 


[7] 


where  the  product  (u  v)  is  calculated  from  each  sounding.  On  the  right 

B  V 

hand  side  of  Eq.  7,  the  first  term,  -  “,  represents  the  total 
effects  due  to  the  mean  radial  circulation,  the  standing  eddy,  and  the 

transient  eddy;  the  second  term,  -  -tallil,  represents  the  effects  due 
to  the  mean  radial  circulation  and  standing  eddy  only,  and  has  been 
included  in  Eq.  5. 

Transient  eddies  have  been  determined  from  Eq.  7.  They  are 
generally  small  at  inner  radii  but  quite  pronounced  at  outer  radii 
especially  at  upper  levels.  For  this  reason,  this  study  will  evaluate 
the  cumulus  convective  effects  only  at  inner  radii  (r  <  5°  latitude) 
where  the  transient  ejjg  effggt  floe?  Ml  Ulitll  ftLflS-t  Ifefi.  results. 
Another  reason  for  limiting  this  study  within  5°  latitude  radius  is  the 
consideration  of  the  generally  observed  radial  extent  of  convection 
associated  with  a  tropical  cyclone  system.  (This  will  be  discussed  in 


the  next  section.) 


Figure*  4*  and  4b  illustrate  the  total  effect.  -  and  the 

mean  aud  standing  eddy  effect,  -  of  1-5°  latitude  radial  and 

axiauthal  averages  for  the  typhoon  (data  set  VPD4)  and  the  hurricane 
(data  set  ATD4).  The  difference  between  the  two  profiles  represents  the 

transient  eddy  effect,  -  This  is  generally  insignificant 

except  at  700  nb  for  1PD4  and  at  upper  levels.  The  transient  eddy 

effect  ten*  in  Eq.  6  (-  Magnitude  of  about  3  ■  s  *day  * 

or  less  at  upper  levels.  Compared  to  the  calculated  residual  CcuF^J 
vhich  has  a  Magnitude  of  about  10  a  s  *day  *  at  upper  levels  (as  will  be 
shown  later)  this  transient  eddy  effect  can,  to  a  first  approximation, 
be  neglected. 


Fig.  4.  1-5°  latitude  radial  averages  of  the  total  effect  of  -  uv  It, 

(solid),  and  aean  and  standing  eddy  effect,  -  [u]  (vl/r. 
(dashed)  for  (a)  typhoon  (VPD4)  and  (b)  hurricane  -[u][v)/r 
(ATD4)  data  sets. 


Although  the  1-5°  averaged  value  of  [u'v']  can  be  quantitatively 
estimated  from  Eq.  7,  the  results  are  not  accurate  enough  to  provide  a 


reliable  estimate  of  the  radial  gradient  term,  ^ •  However,  tak¬ 

ing  the  1-5°  latitude  radial  average  of  this  term  yields 

([u'v']  -  [u'v']  )/4°.  It  is  likely  that  this  term  has  the  same 

5°  1° 

(or  smaller)  magnitude  as  that  of  [u'v']/r  term,  thus  is  probably 


negligible  too.  This  term  is  thus  assumed  to  be  negligibly  small. 


After  neglecting  horizontal  eddy  effects 


the  cumulus  friction  can  then  be  approximated  by 


[cu  F0] 


~  &&  +  SsMhl  +  +  7^}-  +  IsIIvi  _  [W- 

~  dt  xdt  dp  r© 


This  quantity  is  termed  the  residual  cumulus  friction  -  cumulus  friction 
calculated  as  a  residual  from  the  grid  scale  tangential  momentum  budget. 
3.3  Calculated  Residual  -  Cumulus  Friction 

Consideration  of  cumulus  effects  in  a  tropical  cyclone  requires  an 
examination  of  the  distribution  of  cumulus  clouds.  Besides  the  BCE's 
shown  in  Fig.  3,  the  radial  extent  of  total  cloudiness  associated  with 
different  intensity  stages  of  the  tropical  cyclone  in  the  northwest 
Pacific  in  four  cardinal  directions  is  given  in  Table  3  (Arnold,  1977). 
The  high  standard  deviations  indicate  a  substantial  variability  in 
overall  cloudiness.  It  must  also  be  remembered  that  this  mean  radius 
may  enclose  large  cloud-free  regions.  Nevertheless,  these  data  show 
that  most  of  the  convection  associated  with  the  tropical  cyclone  occurs 
at  radii  inside  of  5°  latitude.  Because  very  few  observations  are 


20 


TABLE  3 


Radial  extent 

of  total 

storm 

associated  cloudiness  (after  Arnold,  1977) 

STORM  STAGE 

I 

II 

III 

IV 

Mean  Minimum 
Central 
Pressure  (mb) 

-1005 

-998 

-988 

-959 

Direction: 

Ave . 

SD 

Ave . 

SD 

Ave . 

SD 

Ave . 

SD 

North 

4.6 

2.1 

3.9 

1.5 

6.1 

2.0 

5.3 

2.2 

East 

5.6 

2.7 

5.2 

2.1 

5.9 

2.7 

5.7 

2.5 

South 

5.1 

2.8 

4.7 

2.9 

6.5 

3.1 

5.8 

2.5 

West 

5.0 

2.6 

3.5 

1.8 

5.2 

2.3 

4.9 

1.8 

R 

5.0 

2.6 

4.3 

2.2 

6.0 

2.6 

5.5 

2.3 

Average  (Ave.)  and  Standard  Deviation  (SD)  values  given  in  °Latitnde. 
R  is  the  average  radial  extent  of  four  cardinal  directions. 


available  inside  1°  latitude  radius,  the  residual  cumulus  friction  being 
presented  in  this  paper  represents  the  1-5°  latitude  radial  and 
azimuthal  averages  about  the  cyclone. 

Assuming  that  cumulus  clouds  only  rearrange  horizontal  momentum  in 
the  vertical,  that  is 

lOOmb  _ 

/  Ccu  Fa]dp  =  0, 

P  y 

sfc 

a  correction  factor  must  be  applied  at  each  level  so  that  zero  net 
tropospheric  momentum  change  results.  This  correction  factor  is  just 
the  vertical  mass-weighted  average  of  the  calculated  residual  (see 
section  3.4  for  further  discussion).  The  sign  convention  is  that  a 
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positive  cumulus  friction  will  cause  a  cyclonic  acceleration  (or 
anticyclonic  deceleration)  of  the  wind. 

Figures  5  to  12  illustrate  the  vertical  profiles  of  1-5°  latitude 
radial  and  azimuthal  averages  of  residual  cumulus  friction  for  all  16 
data  sets.  Although  these  profiles  may  be  somewhat  in  error  due  to  the 
absence  of  horizontal  eddies,  uncertainty  in  the  calculation  of  surface 
friction,  noise  inherent  in  the  composite  data,  etc.,  they  do  show 
systematic  variations  between  data  sets  and  are  believed  to  be  reliable 
estimates  of  cumulus  friction.  For  convenience,  1-5°  latitude  radial 
and  azimuthal  averages  will  be  omitted  (but  implied)  in  all  the 
following  discussions,  unless  otherwise  mentioned. 

Figure  5  shows  the  cumulns  friction  profiles  for  different  stages 
of  the  developing  tropical  cyclone  in  the  northwest  Pacific.  These 
profiles  all  have  the  same  pattern:  pronounced  positive  cumulus- induced 
acceleration  at  upper  and  lower  layers,  and  negative  acceleration  at 
middle  layers.  The  positive  acceleration  at  700  mb  may  be  due  to  the 

r'a'w'*  1 

trims lent  eddy  effect,  r  A  (refer  to  Fig.  4).  However,  the 
physical  reason  of  causing  these  spikes  is  still  not  understood.  These 
profiles  also  reveal  that  cumulus  friction  is  directly  related  to  system 
intensity.  This  is  because  the  stronger  weather  systems  have  a  larger 
vertical  wind  shear  (and  probably  more  vertical  mass  recycling)  although 
not  much  difference  exists  in  the  mean  vertical  motion.  The  cumulus 
friction  in  the  lower  layers  for  Atlantic  systems  is  not  as  positive  as 
that  for  the  Pacific  cases  (Fig.  6). 

Figures  7  and  8  illustrate  cumulus  friction  profiles  for  developing 
and  non-developing  cloud  cluster  systems.  The  two  non-developing  cases 
(WPN1,  ATN1)  exhibit  a  similar  pattern:  positive  acceleration  from  the 


pressure  (100  mb) 


Fig.  7. 


The  sane  as  Fig.  S  but  Fig. 
for  developing  (WPD2) 
and  non-developing  (WPN1) 
cloud  clusters  in  the 
northwest  Pacific. 


8.  The  same  as  Fig.  5  but 
for  developing  (ATD1) 
and  non-developing  (A1TJ1) 
cloud  clusters  in 
the  northwest  Atlantic. 


Fig.  9.  The  aaae  as  Fig.  5  but 
for  developing  (ATD2) 
and  non-developing  (ATN3) 
tropical  depressions 
in  the  northwest  Atlantic. 


Fig.  10.  The  same  as  Fig.  5  but 
for  deepening  (WPD3(D)) 
and  non-developing  (WPN3) 
tropical  atoms  in  the 
northwest  Pacific. 
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surface  up  to  600  ab  (slightly  negative  at  900-9S0  ab  for  WPN1),  and 
negative  aloft  except  at  200-250  ab  where  cumulus  friction  is  less 
negative  or  even  positive  (for  VPN1).  The  two  developing  cases  (WPD2, 
ATD1)  also  show  a  siailar  pattern  which  is  qualitatively  different  froa 
the  non-developing  cases.  Although  these  differences  are  within  the 
liaits  of  possible  error,  they  do  show  consistency  for  siailar  stora 
stratification  in  both  ocean  basins. 

One  aight  expect  that  the  difference  between  the  developing  and 
non-developing  cases  would  becoae  larger  as  the  system’s  intensity 
becoaes  stronger.  Figures  9  and  10  show  this  progressive  difference. 
For  the  interaediate  systeas  of  tropical  depression  (stage  2)  or  stora 
(stage  3)  no  apparent  difference  is  found  between  the  intensifying  and 
non-intensifying  cases  at  upper  levels.  Major  differences  do  occur 
below  400  ab,  however.  Cuaulus  friction  tends  to  increase  the  surface 
to  400  ab  vertical  wind  shear  for  the  intensifying  cases  and  decrease 
greatly  this  vertical  wind  shear  for  the  non- intensifying  cases  (ATN3, 
VPN3).  For  the  aature  typhoon  and  hurricane  data  sets  (Figs.  11  and 
12),  cuaulus  friction  increases  the  surface-350  ab  vertical  wind  shear 
nore  for  deepening  cases  (WPD4(D),  ATD4(D))  than  for  filling/steady 
cases  (WPD4,  WPD4(F),  ATD4).  Also,  the  deepening  cases  have  a  smaller 
negative  cuaulus  friction  at  800-900  ab. 

3.4  Discussion 

Cuaulus  friction  has  been  calculated  as  a  residual  of  the 
tangential  momentum  equation  for  16  composite  data  sets.  Remarkable 
consistencies  have  been  found  between  these  profiles:  Pacific  versus 
Atlantic,  strong  versus  weak;  developing  versus  non-developing,  etc. 
Although  these  profiles  are  considered  to  represent  primarily  cuaulus 
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(ms-1  day'1) 


Fig.  11.  1-5°  latitude  radial  and  azimuthal  average*  of  cumulus  fric¬ 

tion,  [cu  F  ] ,  for  mature  (WPD4),  deepening  (WPD4(D))  and 
filling  (WPD?(F))  typhoons  in  the  northwest  Pacific. 


(m»M  day'1) 


Fig.  12.  1-5°  latitude  radial  and  axiauthal  average*  of  cuaulua  fric 

tion,  [cu  F_] ,  for  mature  (ATD4) ,  and  deepening  (ATD4(D) ) 
hurricane*  in  the  northwest  Atlantic. 


friction,  one  has  to  remember  that  they  may  include  other  sub-grid  scale 
effects. 

Recalling  that  all  these  cumulus  friction  profiles  have  been 
adjusted  to  prevent  any  generation  or  destruction  of  tangential  wind, 
cumulus  convection  is  thus  considered  to  be  only  a  momentum 
rearrangement  mechanism.  The  correction  factor,  which  is  subtracted  out 
from  the  directly  calculated  residuals,  represents  the  vertical  mass- 
weighted  average  of  momentum  source  or  sink  if  momentum  generation  or 
destruction  is  present  (a  positive  correction  factor  means  a  generation 
of  cyclonic  momentum).  One  might  speculate  that  such  generation  of 
momentum  is  correlated  to  non-hydrostatic  pressure  gradient  forces 
associated  with  the  activity  of  squall  line  convection  which  somehow 
requires  an  environmental  vertical  wind  shear.  Figure  13  does  show  a 
positive  correlation  between  the  correction  factor  and  the  900-200  mb 
vertical  wind  shear.  However,  this  result  needs  more  verification  in 
other  data  sets  and  more  theoretical  backing  before  we  can  confidently 
conclude  that  cumulus  convection  can  generate  or  destroy  momentum.  At 
this  time  the  author  believes  it  is  more  reasonable  to  treat  cumulus 
convection  as  a  vertical  rearranging  process  only. 


The  relationship  between  the  vertical  wind  shear  n~V_00) 

and  the  correction  factor  (CF  -  vertical  surface  to  100  mb 
mass  weighted  average  of  unadjusted  cumulus  friction),  which 
is  subtracted  out  from  the  calculated  residual  to  assure  that 
the  tropospheric  average  of  cumulus  friction  is  zero. 

Positive  correlation  is  found  between  these  two  parameters. 
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4.  PARAMETERIZATION  OF  CUMULUS  FRICTION 

Although  the  calculated  residuals  are  considered  to  be  reliable 
estimates  of  cumulus  convective  effects,  they  cannot  be  properly 
included  in  the  governing  equations  without  some  parameterization 
scheme.  To  this  end,  many  cloud  models  have  been  developed.  Such  cloud 
models  treat  the  cumulus  clouds  as  an  ensemble  and  vary  the  entrainment 
and  detrainment  rates  for  different  cloud  types.  Residual  cumulus 
friction  calculations  are  used  to  verify  the  cloud  model  diagnoses  (Reed 
and  Johnson,  1974;  Shapiro,  1978;  Cho  t_t  jtl. ,  1979;  Shapiro  and  Stevens, 
1980).  In  Shapiro  and  Stevens  (1980),  three  cloud  models  are  employed 
to  diagnose  the  cumulus  friction;  they  include  spectral,  bulk  and  single 
cloud  models.  They  concluded  that  a  single  cloud  model  can  reasonably 
well  diagnose  the  sub-grid  scale  residuals  (called  apparent  momentum 
source  in  their  paper)  for  tropical  waves  in  the  Atlantic  GATE  region. 

In  this  study,  a  simple  single  cloud  model  is  developed  to  diagnose 
cumulus  friction  in  the  region  1-5°  latitude  radius  about  the  mature 
typhoon  and  hurricane.  The  model-diagnosed  results  are  then  compared  to 
the  residual  cumulus  friction  presented  in  Chapter  3.  Although  a 
consistency  between  the  calculated  residual  and  the  model-diagnosed 
cumulus  friction  may  not  provide  a  better  understanding  of  mesoscale 
dynamics,  it  does  give  certain  confidence  as  to  how  cumulus  friction 
effects  can  be  parameterized  in  numerical  models.  Because  of  our  lack 
of  knowledge  of  mesoscale  dynamics,  it  is  not  necessarily  true  that  a 
more  complicated  cloud  model  can  handle  this  cumulus  friction  better 
than  a  simple  model  can.  In  addition,  any  consistency  between  the 
observed  residual  and  the  model  results  can  be  considered  as  general 


support  to  the  hypothesis  that  cuotluj  convective  transports  can  be 
calculated  as  a  residual. 

4.1  Formulation 

Total  convective  effect  (denoted  as  TCP.)  is  defined  as  including 
both  verticil  near!  and  eddy  transports.  This  can  be  decomposed  into 
upward  and  downward  motion  effects  as: 

TCE  =  (-  +  (- 

dp  dp 


3u>  v 

=  (-  -?-*>  +  c  -  — M> 


3p 


3p 


[9] 


where  subscripts  u  and  d  denote  upward  and  downward  motion  respectively. 
Because  the  tropical  atmosphere  is  so  stable,  all  the  upward  motions 
have  to  bs  associated  with  cloud  actions.  For  simplicity,  a  single 
cloud  model  night  be  utilized  where  the  do»r.ward  motion  is  considered  to 
be  the  gentle  broadacale  subsidence  having  the  same  momentum  as  the 


environment.  Equation  9  then  becomes 


TCE  =  - 


3w  v  3u>  ,v 
c  c  _  _ d  e 

3p  3p 


[10] 


where  ei  and  <o.  are  the  average  cloud  and  downward  mass  fluxes  over  a 
c  d 

grid  spacing  respectively;  vc  and  are  the  cloud  and  environmental 

momentuns  respectively.  Values  of  u>  ,  v  ,  u>,  and  v  are  taken  to  be 

cede 

constant  at  any  level.  Since  [u]  =  u>  + 

c  d 


0).  =  M  -  u>  . 

o  c 


[11] 


Also, 


[▼]  =  «tv  +  (l-a)v 

c  e 


[12] 


where  a  is  the  fractional  area  covered  by  cumulus  clouds.  If  cumnlns 


clouds  only  occupy  •  very  small  fraction  of  the  total  area,  i.e.»  o  << 


1,  Eq.  12  becomes 


v  -  [▼]. 

© 


Substituting  (11)  and  (13)  into  (10)  gives 


Thus, 


"  Tp  KVc  +  ([wl  "  \)ly]\ 


mean  convective  eddy  convective 

effect  effect 


and  cumulus  friction  can  be  parameterized  as 


r  c  i  _  _  v  1 

lcu  FG]  *  dp 

Or,  talcing  the  azimuthal  average  yields 


9w  (v  ~[vl) 
c  c _ 


[cu  Fel  -  - 


9u>  (v  ~[v]  ) 
c  c _ 


It  may,  however,  be  inaccurate  to  assume  that  all  the  downward 

motions  are  in  the  form  of  gentle  broadscale  subsidence.  Part  of  the 

downward  motion  might  be  in  the  form  of  strong  cloud  downdrafts  or 

mesoscale  downdrafts  (Johnson,  1976,  1980;  Zipser,  1977)  which  could 

conceivably  carry  upper  level  momentum  down  to  the  lower  level.  Then 

Eq.  9  should  likely  be  rewritten  as 

9w  v  9ii>.  v.  die.  v. 

ICE  =  -  -  Z  U6] 

dp  dp  op 

where  the  subscript  dg  denotes  the  gentle  subsidence,  and  dn  any  other 
kind  of  strong  downdrafts.  Through  the  same  argument  as  before,  Eq.  16 


becomes 


[17] 


,  Ou)  (v  -[%•])  <:  j!>  (J  .  (V  -  [  v] )  I 

_  _  0  I  u< )  t  v  j  _  c  c _ _  _  t  p  a  _ 

dp  bp  bp 


end  the  azimuthal  averaged  cumulus  friction  beccaes 


Ccu  V  =  '  "  c 


OUJ  (  V  -  [vj  ) 
c  c 


?!l>  X-.<V.  -  i  V  J )  [ 
_<  u  d  du _ 

dp 


[18] 


where  D  -  (flu,  j/[ie.])  is  the  fraction  of  local  dowr.'-ard  mass  flux 
n  dn  a 

which  is  accomplished  by  the  nth  kind  of  s-iv/Cg  downJraft.  Equation  18 
is  identical  to  Eq.  15  if  there  is  no  otner  kind  of  string  downdraft 

present,  i.e,,  D  =0.  Note  that  in  both  Eqs.  15  and  18  the  momentum 

n 

carried  Iowa  by  gentle  brosdscale  subsidence  is  assumed  to  be  the  same 

as  the  grid  box  average  momentum  [\ ]  (rq.  12).  fhis  assumption  might 

not  hold  well  at  lower  levels  where  the  cumulus  clouds  ‘night  occupy  up 

to  15-20‘iu  of  the  total  area  such  that  o  is  not  much  less  than  one. 

4.2  .Moisture  Budget  and  Vertical  Mass  Recycling 

Hie  determination  of  vertical  mass  recycling  (or  a>  )  is  required  to 

c 

diagnose  cumulus  friction  from  Eqs.  15  and  18.  Frank  (1977c)  employed  a 
spectral  cloud  model  and  used  the  moist  static  energy  (h)  budget  (first 
discussed  by  Tanai  ejt  .al.  (1973)  and  Lopez  (1973))  to  diagnose  the 
vertical  mass  recycling,  but  strong  downdrafts  were  not  considered  in 
these  studies.  Inclusion  of  cloud  downdrafts  or  mc3osoale  downdrafts 
can  cause  important  changes  in  vertical  mass  recycling  (Johnson,  1976, 
1980).  To  a  large  extent,  the  type  of  model  to  be  employed  determines 
the  magnitude  of  the  vertical  mass  recycling  one  obtains.  In  this  stndy 
a  single  cloud  model  is  employed,  rather  than  a  more  complex  spectral 
cloud  model.  This  enables  s  simpler  physical  interpretation  of  the 


results 


Following  Gray  (1973)  for  a  steady  state  cloud  cluster  system,  the 
required  upward  moisture  flux  carried  by  the  sub-grid  scale  circulation 
at  any  level  k  (k  =  0  at  surface)  is  given  by: 


lE.t.  "  'Vi-il  -  IS  »'  (C-E)kl 


Required  upward 

/Surface  evaporation 

Condensation  to  rain 

aoisturo  flux 

1 (k=l ) 

or  Q  at 

or  condensation  minus 

carried  by 
local  circula¬ 
tion  at  level  k 

{level 

k-1  fk  >  1) 

reevaporation  between 
level  k  and  k-1. 

+  q)k  -  (<*>  qt> J  +  {-  (\7-Vq)^  ^  | 


[19] 


Vertical  moisture 
convergence  carried 
by  mean  vertical 
motion  between 
level  k  and  k-1. 


Horizontal  moistnre 
convergence  between 
level  k  and  k-1. 


where  k  denotes  the  layer  between  level  k  and  k-1,  and  the  overbar 
(  ),  denotes  the  mean  circulation  (or  azimuthal  average).  Horizontal 
transport  by  transient  eddies  is  neglected.  The  major  problems  in  this 
calculation  are  the  determination  of  surface  evaporation  and  the  amount 
of  condensation  that  goes  to  rain  at  each  layer. 

I 

4.2.1  Surface  Evaporation 

As  discussed  by  Gray  (1981)  and  Fingerhut  (1982),  use  of  the  bulk- 
formula  to  calculate  surface  evaporation  in  the  high  wind  speed 
condition  of  the  tropical  cyclone  may  result  in  a  significant 


under e  s t i ma  t ion. 


If  reasonably  accurate  rainfall  data  could  be 
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obtained,  it  could  be  used  to  calculate  surface  evaporation  as  a 


residual  from  the  formula 


(Surface 
evapor a  t ion 


(Total  horizontal  moisture 
convergence  from  surface 
to  tropopause 


[2 


where  no  net  moisture  change  (in  time)  is  assumed  through  the  vertical 
column.  Using  rainfall  observations  from  16  hurricanes  over  Florida, 
Miller  (1958)  calculated  rainfall  amounts  associated  with  hurricanes. 
Topographic  effects  were  assumed  to  be  negligible  in  these  data  due  to 
the  flat  terrain  of  Florida.  Measurement  errors  caused  by  high  winds 
were  also  corrected.  The  rainfall  amounts  associated  with  typhoons  hav 
been  obtained  by  Frank  (1976)  using  the  observed  rainfall  data  from  13 
small  islands  in  the  northwest  Pacific.  Figure  14  illustrates  these 
rainfall  amounts  associated  with  hurricanes  as  well  as  typhoons.  Very 
reasonable  consistency  is  found  between  these  two  curves. 

Using  these  rainfall  data,  the  surface  evaporations  of  mature 
hurricanes  and  typhoons  at  1-5°  latitude  radius  about  the  system  centei 
are  determined.  Results  are  shown  in  Table  4.  They  are  quite  similar 
to  those  obtaued  by  Frank  (1976)  for  the  typhoon.  These  results  also 
reveal  that  a  typhoon  system  has  almost  the  same  surface  evaporation  a 
that  of  a  hurricane  system  at  1-5°  latitude  radiu*.  The  greater 
rainfall  amount  associated  with  the  typhoon  system  is  •  result  of  larg 


horizontal  moisture  convergence. 
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Fig.  14.  Observed  precipitation  around  tropical  cyclones  from  Florida 
and  from  N.W.  Pacific  Islands  (after  Frank,  1976). 


TABLE  4 

Vertical  integrated  moisture  budget  associated  with  the  typhoon  and 
hurricane  at  1-5  latitude  radius  about  the  system  center  in  cm/day. 


Typhoon 


Hurricane 


Observed  Rainfall  2.9  2.5 

Horizontal  Moisture  Convergence  1.4  0.9 

Required  Surface  Evaporation _ 1.5 _ 1.6  _ 

4.2.2  Vertical  Distribution  of  Condensation  to  Rain 

An  ad-hoc  approach  was  used  to  determine  the  condensation  to  rain, 

or  condensation  minus  re-evaporation  (C-E) ,  denoted  as  R  at  each  layer. 

The  rainfall  R  is  decomposed  into  two  parts,  one  part  from  mean  vertical 

motion  (R^)  and  the  other  part  from  local  recycling  (R^).  The  attempt 

here  is  to  make  the  ratio  of  the  local  recycling  (»  )  to  the  mean 

vertical  motion  close  to  the  ratio  of  recycling  rainfall  (Rf)  to  mean 
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I 

I 


vertical  motion  rainfall  (R  )  8t  each  layer  or 

m  * 


(u  /w)  ~  (R  ,/R  ) 

r  k  ~  r  m  k 

In  addition*  the  resulting  vertical  cass  fluxes  have  to  resemble  those 

obtained  by  Frank  (1977c)  when  no  strong  Jov ndrsfts  are  considered. 

The  condensation  due  to  rre&n  vertical  action  (C  )  can  be  calculated 

m 

as ; 

3  u>  q  , 

C  =  -  (— r - *)  ^ 

a  3p  g 

Assuming  half  of  the  liquid  water  condensed  due  to  mean  upward 
circulation  is  re-evaporated  (E  )  then 


or 


R  =  C  -  E  =  0.5  C 

m  an  d 


(Vk  =  Y~g  !(w  Vk  -  (U  qs)k-l| 


The  total  condensation  to  rain  due  to  the  local  recycling  ((Z  R  ).  )  can 

k 

then  be  obtained  as  a  residual  by 


Z  (R  ).\ 
t  ' ») 


-  (Total  rainfall)  -  /Z(R  >. 

(k  “  1 


The  vertical  distribution  of  Rf  i*  determined  by  the  mixing  ratio 

(q)  and  the  difference  between  the  saturation  mixing  ratio  and  the 

mixing  ratio  (Aq  =  qs~q) »  or 

j  |  (qAp)^  (AqAp)^ 


(Vk  "  2  )  Z(qAp)“  +  Z(AqAp) 

v  k  -  k 

Therefore,  the  condensation  to  rain  at  each  layer  is 


k  M} 


(R  ).  +  (R  >.. 

m  k  r  k 


[21] 


The  vertical  distribution  of  R,  R_  snd  R  sre  shown  in  Table  5. 

D  X 


TABLE  5 


Condensation  to  rain  (mm/dey)  due  to  mean  vertical  motion  (R  ) ,  local 
recycling  (R  )  and  total  upward  motion  (R)  for  the  typhoon  ($PD4)  and 
the  hurricane  (ATD4)  at  1-5  latitnde  radius  about  the  system  center. 

Typhoon  (tfPD4)  Hurricane  (ATD4) 


Layer (mb) 

R 

m 

R 

r 

R 

R 

m 

R 

r 

R 

100-200 

0.07 

0.03 

0.1 

0.03 

0.02 

0.05 

200-300 

0.73 

.33 

1.06 

0.46 

0.33 

0.79 

300-40C 

1.47 

1.03 

2.50 

1.13 

0.77 

1.90 

400-500 

1.51 

1.32 

3.33 

1.22 

1.40 

2.62 

500-600 

1.40 

2.55 

3.95 

1.17 

2.02 

3.19 

600-700 

1.49 

3.36 

4.85 

1.24 

2.72 

3.96 

700-800 

1.23 

4.00 

5.23 

1.09 

3.43 

4.52 

800-900 

0.96 

4.40 

5.36 

0.97 

4.06 

5.03 

900-950 

0.46 

2.51 

2.97 

0.43 

2.38 

2.81 

Total 

9.32 

20.03 

29.35 

7.74 

17.13 

24.87 

4.2.3  Requi 

red  Local 

Vertical 

Moisture  Flux  and  Mass 

Recycling 

Equation  19  then 

can  be  fully  evaluated  to  determine  the  required 

vertical  moisture  flux 

.  This 

is  shown  in 

Fig.  15  for 

the  typhoon 

and 

the  hurricane.  It  can 

be  seen 

that  these 

two  curves 

are  almost 

identical.  The  required  vertical  mass  recycling  («  )  can  be 

determined  from  (Gray,  1973) 

<“r>k  -  -  («r)k  Vk  -  law 

where  subscript  u  and  d  denote  upward  and  downward  as  before.  The 
upward  cloud  mass  flux  is  given  by: 
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Qf  (mrn  day'1) 


Fig.  15.  Required  vertical  moisture  flux  carried  by  local  circulation 
(Qf)  for  it‘PD4  (solid)  and  .'.7D4  (cashed)  datB  sets, 

u>  =  tn  +  e>  [23] 

c  r 

and  the  downward  mass  flux  by: 

us,  -  -tn  [24] 

c  r 

The  quantity  q  .  can  be  approximated  by  q  (calculated  from  mean 

temperature)  if  all  the  upward  motions  are  assumed  to  be  saturated. 

This  approximation  is  affected  by  two  factors:  1)  sene  upward  motion 

might  be  sub-saturated  and  2)  a  slightly  higher  temperature  might  be 

expected  inside  the  convective  clouds.  However,  these  two  effects  will 

tend  to  cancel  each  other.  If  the  downward  motion  is  assumed  to  be  in 

fora  of  gentle  subsidence,  then  by  the  same  argument  as  that  of  the 

momentum,  q.  can  he  approximated  by  q  •  The  profile  of  (q  -q)  is 
o  c  * 

shown  by  the  solid  line  in  Fig.  16.  The  required  local  mess  recycling 
was  then  calculated  (see  Fig.  17).  The  close  similarity  between  the 


Fig.  16.  Assumed  mixing  ratio  anogalies  for  cloud  updraft*  (q  -  q) , 


cloud  downdrafts 
(%,d  "  qK 


(Q  *  ~  q)  and  meso-scale  downdrafts 
cd 
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Fig.  17, 


Mean  vertical  mas*  flux  (w),  cloud  mass  flux  (»,)  and  down 
ward  mass  flux  (u  )  for  WPD4  (solid)  and  ATD4  (dashed)  data 
sets  when  no  strong  downdraft  is  included. 


tyhooon  and  the  hurricane  is  again  seen.  Ti. •.  :■  c  .  r;:  i'.-s  are  very 
similar  to  those  obtained  by  Frank  (i?77c). 

4.2.4  The  Effect  of  Downdrafts  on  the  Diagnosed  Local  Mass  Recycling 

As  mentioned  before,  a  part  of  the  downward  mass  flux  night  be  in 
the  form  of  strong  downdrafts  which  have  considerably  different 
thermodynamic  aspects  than  these  of  large  scale  gentle  subsidence. 
Therefore,  two  kinds  of  downdrafts  other  then  gentle  subsidence  are 
considered,  following  Johnson  (le"6,  .JhO):  (a)  convective  scale 
downdrafts  from  350  r.b  (entrainment  rate  X  ~  0.05  be  ^  in  Johnson's 

model,  1976);  (b)  nesoscale  downdrafts  from  600  r.’o.  The  difference  in 
mixing  ratio  between  the  convective  downdrafts  (values  .-.re  taken  from 
Johnson,  1976)  and  the  c  nv  iron-vent  ( approx  i~.a  ted  by  q)  for  the  1-5° 
latitude  radius  average  of  a  typhoon  is  shv.1.  by  the  decked  line  in  Fig. 
16.  The  profiles  for  t^e  hurricane  are  approximately  the  same  as  those 
for  the  typhoon.  This  convective  downdraft  is  cool  ar.d  dry  when  it 
reaches  the  surface  due  to  re-evapora tion.  An  increase  of  vertical  mass 
recycling  is  expected  if  this  convect ive-scr 1 e  downdraft  is  considered 
because  it  is  more  moist  than  the  environmental  air. 

Table  3  of  Johnson's  (1980)  paper  shows  the  average  thermodynamic 
structure  of  mesoscile  dc  -virafts  obtained  from  Zipser's  (1977) 
composite  standings.  The  mixing  ratio  for  this  ciesoscale  downdraft  was 
reproduced  from  this  table.  The  dotted  line  in  Fig.  16  shows  the  mixing 
ratio  difference  between  the  aesoscale  downdraft  and  the  environment  for 
the  typhoon  data  set  (WPD4).  This  profile  shows  that  mesoscala 
downdrafts  are  much  drier  than  the  environment,  and  would  thus  act  to 
reduce  the  amount  of  required  vertical  mass  recycling  for  water  vapor 


balance 


0 


For  simplicity  (also  due  to  the  lack  of  good  quantitative  knowledge 
concerning  these  downdrafts),  the  average  mixing  ratio  associated  with 
all  the  downward  motion  i«  teken  to  he:  (a)  q  ,  if  p  <  300  mb; 

(b)  (  q  +  qc(j)/2,  if  300  nb  <  P  i  600  +  9cd  +  Sod^3,  i{  P 

>  600  mb;  where  subscripts  cd  and  md  represent  cloud  downdraft  and 
aesoscale  downdraft  respectively.  The  vertical  r-.ass  recycling  at  each 
level  is  thus  given  by 

(id  )  -  -  Q  g/(q  -  q  )  [25] 

x  r  g  a 

and  cloud  and  dcvnward  mass  fluxes  are  determined  by  Eqs .  23  and  24. 

These  results  are  shown  in  Fig.  18  where  the  low  level  vertical  mass 

flux  is  slightly  reduced.  It  must  be  remembered  that  the  profiles  shown 

in  Fig.  18  will  certainly  vary  from  reality  due  to  the  assumptions 

involved  in  this  simple  model.  However,  these  profiles  are  still 

physically  reasonable  and  show  consistency  for  both  ocean  basins. 

As  mentioned  before,  the  ratio  of  the  mean  vertical  motion  to  the 

calculated  local  recycling  has  to  match  the  ratio  of  the  assumed 

condensation  to  rain  produced  by  the  mean  vertical  motion  and  the 

condensation  to  rain  produced  by  the  local  recycling.  Figure  19a 

depicts  the  assumed  profiles  of  condensation  to  rain  for  the  three 

quantities  of  R  ,  R  and  R  at  each  level.  Note  that  these  profiles 
m  r 

resemble  those  of  <e,  o>  and  u>  quite  well.  It  should  be  noted  that 

r  c 

the  calculated  local  recycling  profile  is  very  sensitive  to  the  assumed 

vertical  distribution  of  R  (or  R) .  If  a  larger  R  is  assumed  at  lower 

r  r 

levels,  the  magnitude  of  w is  reduced  substantially  at  all  the  higher 
levels.  For  comparison  purposes,  the  ratios  of  the  assumed  R  to  R  and 

IB 

of  the  id  to  the  calculated  w  are  shown  in  Fig.  19b.  These  two 

c 

profiles  are  almist  identical  at  middle  to  lower  le\cls.  At  upper 


pressure  (lOO  mb) 


( loo  mb) 


Fig.  18.  Mean  vertical  mass  flux  (u>)»  cloud  mass  flux  (u>  )  and  down¬ 
ward  mass  flux  (oi^)  for  «iJD4  (solid)  and  ATD4  (dashed)  data 
sets  when  both  ki:ds  of  downdrafts  are  included. 


(a)  (b) 

Fig.  19.  (a)  Condensation  to  rain  for  every  100  mb  layer  due  to  mean 

vertical  mass  flux  (R  ),  local  mass  recycling  (R  )  and  the 
total  upward  mass  flux  (R).  (b)  Ratios  of  R  to  R  compared 

to  that  of  mean  vertical  mass  flux  (  w)  to  tSe  total  upward 
mass  flux  (u>  ) . 


levels,  the  difference  is  ft  little  larger.  Th  i s  is  because  the 
condensation  magnitudes  are  much  smaller  at  upper  levels  and  a  small 
calculation  error  can  cause  larger  percentage  differences.  The 
agreement  of  these  profiles  provides  good  support  to  this  mass  recycling 
model . 

4.3  Tangential  Momentum  Associated  With  Upward  and  Dov.n-.tard  Motions 

Besides  upward  and  downward  mass  fluxes,  the  tangential  momentum 
associated  with  cloud  uj'drafts  and/or  downdrafts  is  still  required  in 
order  to  fully  evaluate  Eq.  IS  or  Eq.  18.  It  is  likely  that  the  cloud 
momentum  should  vary  as  the  ascending  air  parcel  entrains  and/or  mixes 
with  the  environ. ental  air.  Since  deep  cumulus  convection  has  a  very 
short  parcel  t  1  ne  scale,  this  vertical  variation  of  cloud  momentum 
should  be  negligibly  small  (Schneider  and  Lindzen,  1976;  Stevens  et  al . . 
1977).  Indeed,  the  diagnosed  cloud  momentum  of  Shapiro  and  Stevens 
(1980)  appears  to  he  fairly  constant  with  height  with  a  magnitude  close 
to  the  environmental  momentum  at  cloud  base  height  (950  mb).  The  cloud 
momentum  values  assumed  for  this  study  are  shown  by  the  solid  lines  in 
Figs.  20a  and  20b.  These  cloud  momentum  values  are  assumed  to  vary  in 
the  vertical  as  the  ascending  air  entrains  environmental  air.  A  small 
constant  mixing  rate  (both  entrainment  and  detrainment)  between  clouds 
and  the  env  i  i  .-i-iH-nt  is  also  used. 

In  the  diagnosis  of  vertical  mass  recycling,  the  effects  of  cloud 
d ovr drafts  and  l-.esoscale  downdrafts  are  also  considered  as  those  in  the 
cumulus  parameterization  scheme  (Eq.  18).  The  vertical  profiles  of 
momentum  associated  with  these  two  kinds  of  downdrafts  are  shown  by  the 
dotted  lines  in  Fig.  20.  The  net  entrainment  and  mixing  effects  again 
are  considered  to  allow  these  downdrafts  to  exhibit  different  momentum 


Fig.  20.  Mean  tangential  wind  profile  ([  v] )  and  tangential  momentum 
associated  with  cloud  updrafts  (V  ),  cloud  downdrafts  (V  ) 
and  meso-scale  downdrafts  (V  )  for  (a)  typhoon  (YPD4)  ana  (b) 
hurricane  (ATD4 )  data  set*.  °d 

characteristics  as  they  descend.  As  discussed  before,  the  momentum 
associated  with  the  broadscale  subsidence  is  considered  to  be  the  same 
as  the  grid-box  average  momentum  {Eq.  13).  Such  an  assumption  is 
somewhat  violated  especially  at  the  lower  levels  where  deep  and  shallow 
convective  clouds  cover  a  significant  percentage  of  the  total  area. 

Although  convective  clouds  are  considered  to  transport 
environmental  momentum  up  and  down,  it  Is  unlikely  that  such  assumed 
momentum  profiles  can  be  observed  locally  in  individual  cases. 
Complications  arise  due  to  the  usuelly  observed  mesoscsle  high  pressure 
associated  with  convective  cells  (Zipser,  1977).  Such  convective- 
induced  pressure  gradients  might  cause  t  substantial  alteration  of  the 


momentum  fields  locally.  However,  it  is  possible  thet  such 


assumed 


momentum  profiles  can  be  obtained  by  compositing  many  observed  momentum 
profiles  associated  with  many  convective  systems.  In  this  composite 
situation,  such  pressnre  gradient  forces  will  probably  average  to  zero. 
4.4  Results  and  Oiscnssions 

The  dashed  lines  in  Fig.  21  show  the  mode l~d i aguosed  cumulus 
friction  without  strong  downdrafts  (Eq.  IS).  The  solid  lines  represent 
the  residual  cumulus  friction.  Remarkable  consistency  was  found  between 
the  model  diagnosis  and  the  calculated  residuals  for  both  the  typhoon 
and  the  hurricane.  The  differences  between  these  two  profiles  for  both 
oceans  are  within  the  limits  of  calculation  error  except  at  700  mb  for 
the  typhoon  case  where  the  difference  appears  too  large  to  be  accounted 
for  by  computation  noise.  This  may  be  due  to  the  fact  that  horizontal 
eddy  effects  were  neglected  in  the  residual  calculations  (see  Fig.  4). 
These  results  are  more  consistent  than  those  (momentum  fields)  obtained 
by  Shapiro  and  Stevens  (1980).  This  is  probably  because  an  evaluation 
of  the  pressure  gradient  force  was  required  in  their  calculations  but 
not  in  this  analysis.  This  may  be  especially  relevant  at  the  upper 
levels  where  pressure  gradient  measurements  are  most  unreliable. 

The  use  of  D  =  1/3  in  Eq.  18  for  both  cloud  downdraft  and 
n 

mesoscale  downdraft  only  changes  the  low  level  results  slightly,  as 
shown  by  the  dotted  lines  in  Fig.  21.  Thus,  although  strong  downdrafts 
have  to  be  considered  in  diagnosing  the  correct  vertical  mass  recycling, 
such  strong  downdrafts  can  be  largely  neglected  in  cumulus  momentum 
parameterization  schemes.  However,  such  a  conclusion  is  somewhat 
dependent  on  the  middle  to  low  level  vertical  wind  shear  which  is  very 
small  in  these  cases.  It  is  possible  that  strong  downdrafts  can  play  a 
much  more  important  role  if  the  middle  to  low  level  wind  shear  is  large. 


Fig.  21.  Model  diagnosed  cumulus  friction  with  (dotted)  and  withont 
(dashed)  both  kinds  of  downdrafts  compared  to  the  residual 
cumulus  friction  (solid)  for  (a)  typhoon  (WPD4)  and  (b) 
hurricane  (ATD4)  data  sets. 


Although  only  a  simple  single  cloud  model  is  used  to  parameterize 
cumulus  friction,  the  model-diagnosed  results  agree  reasonably  well  with 
the  calculated  residuals  in  both  pattern  and  magnitude.  This  is  as  good 
as  one  might  expect  given  the  limitations  of  both  the  observational  data 
and  the  theoretical  model  assumptions.  This  simple  scheme  (Eq.  15)  has 
also  been  used  by  Stevens  e_t  a_l.  (1977)  and  Fingerhut  (1982)  in 
numerical  model  studies  which  yielded  realistic  results. 

Note  that  all  the  discussions  here  are  on  the  net  effect  of  cnmulns 
convection  on  the  large-scale  momentum  fields.  This  net  effect  appears 


to  be  not  strongly  dependent  on  the  detailed  internal  dynamics  of  the 


convective  systems  (also  discussed  by  Cho  e_t  a  1 . ,  1979;  Shapiro  and 
Stevens,  1980).  However,  the  clond  internal  dynamics  must  be  taken  into 
account  shonld  a  single  convective  cell  or  mesoscale  system  be 
considered. 

Shapiro  and  Stevens  (1980)  also  emphasized  the  importance  of  the 
correct  determination  of  the  clond  mass  flux.  However,  although  the 
cloud  mass  flux  diagnosed  in  section  4.2  may  not  be  exact,  realistic 
results  were  still  obtained.  It  is  probable  that  this  parameterization 
scheme  is  more  sensitive  to  cloud  momentum  than  to  cloud  mass  flux. 

More  observational  studies  of  a  spectrum  of  cloud  monentums  are  needed 
to  confirm  fully  these  model  results.  Aside  from  all  these  assumptions 
and  complications,  we  were  nevertheless  gratified  to  learn  that  there 
was  general  consistency  between  the  calculated  residuals  and  the  model 
diagnosed  results.  This  provides  a  good  justification  to  the  hypothesis 
that  the  calculated  residual  at  inner  1-5°  latitude  radii  of  a  tropical 
cyclone  system  is  a  realistic  estimate  of  cumulus  friction. 


Tilt  HOLE  OF  CUMULUS  .MOMENTUM  TRAN  sl*0  RTS 


McBride  and  Zohr  (1981)  found  that  the  most  distinct  difference 
between  developing  and  non-developing  cloud  clusters  is  the 
enviio'  vertical  wind  shear:  the  developing  system  is  embedded  in 

an  environ  tent  v  ith  puch  larger  vertical  wind  shear.  A  zero  vertical 
vir.d  shear  li;:e  extends  across  the  center  of  the  developing  system.  As 
the  system  intensifies  from  a  cloud  cluster  to  a  mature  tropical  storm, 
the  vertical  wind  shear  also  increases  over  a  very  large  radial  domain. 
This  increase  in  vertical  wind  shear  is  largely  due  to  the  increase  of 
mid-low  level  cyclonic  wind  at  inner  radii  and  to  the  increase  of  upper 
level  anticyclone  vind  at  outer  radii  (Fig.  22). 

As  have  been  discussed  by  Silva  Dias  and  Schubert  (1979),  and 
Schubert  et  si.  (1930),  geostrophic  or  cyclostrophic  adjustment  on  the 
cyclone  scale  in  the  tropics  is  primarily  an  adjustment  of  the  mass 
field  to  the  wind  field.  Thus,  for  genesis  and  intensification  to  occur 
a  tropical  system  must  form  an  anticyclonic  vertical  wind  shear  around 

Q 

its  center  (1-5  latitude  radius  region  or  even  further  outward)  to 
sustain  an  upper  level  warm  core.  This  upper  level  warm  core  is 
primarily  attributed  to  the  release  of  latent  heat  (in  form  of 
subsidence  warning)  of  cumulus  convection.  The  understanding  of  how 
this  vertical  wind  shear  can  increase  during  the  evolution  of  the  system 
becomes  fundamentally  important  to  the  tropical  cyclone  genesis  and 


intensification  problem 


Fig,  22.  Vertical  profiles  of  the  mean  tangential  wind  at  (a)  1-5° 
latitude  radial  belt  and  (b)  10  latitude  radius  for  cloud 
cluster  (I),  tropical  depression  or  storm  (II),  and  mature 
typhoon  or  hurricane  (III). 


v  :* 


effect  (ICE) 
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i  di-c  issioas  suggest  that  the  net  effect  of  the 
„•  _ .  -  :  .■!  ci.i'cre  gradient  is  not  important  vhen  averaged  over  a 

- :  t  :..j  -  r  .  J  the  system  center  at  any  radius).  It  is  still 
v.  ,  i  io  examine  the  influences  of  all  the  vertical  motion 

on  the  •  ...  ; ..  ;  eiiium  field  before  the  dovn-grad  lent  or  up-grad i  nt 

■.tit  arc  of  .  it*:  •.  .us  convection  processes  can  be  determined.  In 
o’-  ion  4.1  local  convective  effect  (TCE)  is  decomposed  into  the 

-  -  n  (-  3 1  tr  '  *•  «'t!p)  end  eddy  (-  3[u»*v']/3p)  vertical  transports  (fq. 

.4).  The  r  .  -r.  v.  rticsl  trensp/Ort  can  be  directly  calculated  and  the 
eddy  Ov.  nv.-c  rive  transport  is  approximated  by  the  calculated  residual  in 
the  TCE  otlcuis.  .on. 

The  TCI  profiles  for  the  16  composited  tropical  cyclone  data  sets 
within,  the  1-5"  latitude  radius  region  around  the  system  center  are 
depicted  in  Pips.  23-30.  These  profiles  show  systematic  variations 
similar  to  those  of  the  residual  cumulus  friction,  because  the  mean 
convective  transport  has  the  same  pattern  (negative  acceleration  at 
lower  levels  aud  positive  acceleration  at  upper  levels)  for  all  the  data 
sets.  The  triangle  on  each  curve  shows  the  level  of  maximum  cyclonic 
wind.  It  is  here  where  a  pronounced  negative  acceleration  influence 
(which  acts  to  .pin  down  the  cyclonic  momentum)  is  found  for  every  data 
set.  For  most  of  the  data  sets  this  level  is  aiound  950-800  mb  where 
the  strongest  convergence  is  typically  found.  The  exceptions  are  the 
ton-developing  cluster  data  sets  (WPN1,  ATN1)  which  have  their  maximum 
cyclonic  vinl  at  a  higher  level  (600  mb  and  500  mb).  A  pronounced 
positive  value  (~ii-ir,;um  negative  for  ATN1)  is  found  near  150-300  mb 
maxi  n  divergent  layer  (or  deep  dead  detrainment  layer) 


here  the 
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Fig.  23.  1-5°  latitude  radial  and  azimuthal  averages  of  total  convec¬ 

tive  effect  (TCE)  for  the  northwest  Pacific  (WP)  developing 
tropical  cyclone  systems:  a)  WPD2  (cloud  cluster),  b)  WPD3 
(intensifying  tropical  storm),  c)  WP1)4  (mature  typhoon)  and  d) 
WPD5  (supertyphoon).  The  triangle  shows  the  level  of  maximum 
cyclonic  wind. 


Pig.  24.  1-5°  latitude  radial  and  azimuthal  averages  of  total  convec¬ 

tive  effect  (TCE)  for  the  northwest  Atlantic  developing  tropi¬ 
cal  cyclone  systems:  a)  ATD1  (cloud  cluster),  b)  ATD2  (tropi¬ 
cal  depression),  c)  ATD3  (intensifying  tropical  storm)  and  d) 
ATD4  (mature  hurricane).  The  triangle  shows  the  level  of  max¬ 
imum  cyclonic  wind. 
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TCE  ;rn  s*1  day "') 

The  .  ;?r,s  as  Fig.  23  bat 
for  developing  (V?D2)  and 
non-developing  (WPN1)  cloud 
clusters  in  the  northwest 
Pacific. 
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TCE  i'rri s'1  day*') 

Fig.  26.  The  se.ne  as  Fig.  23  but 

for  developing  ( ATD1 )  and 
non-developing  (A"tNl)  cl  nd 
clusters  in  the  northwest 
Atlantic. 
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WPD3  (D)  i\  \WPN3 
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TCE  jib'1  day”') 

The  same  a*  Fig.  23  but 
for  developing  (ATD2)  and 
r.c  a -dwe  Ic.ping  (ATN 3 ) 
t.ropicH  1  depressions  in  the 
r.  z  r  t  e :  t  Atlantic. 


I  4  8 

TCE  (m  day'1) 

Fig.  28.  The  sane  as  Fig.  23  but 
for  deepening  (WPD3(D)) 
and  non-de vc 1 cp i ng  (  o\ 3 ) 
tropical  stones  in 
the  northwest  Pseific. 
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Fig.  2°.  1-5°  latitude  radial  end  azimuthal  averages  of  total  connec¬ 

tive  effect  (TCE)  for  mature  (WPD4),  deepening  (FPD4(D))  and 
filling  (WI*n4(F) )  typhoons  in  the  northwest  Pacific.  The  tri¬ 
angle  shows  the  level  of  maximum  cyclonic  wind. 


Fig.  30.  1-5°  latitude  radial  and  azimuthal  average*  of  total  convec¬ 

tive  effect  (TCE)  for  mature  (ATD4)  and  deepening  (ATP4(d)) 
hurricanes  in  the  northwest  Atlantic.  Th«  triangle  shows  the 
level  of  maximum  cyclonic  wind. 
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'  icSullS  lead  to  the  CCu-.  2  Us  iou  that  uuuvcglivt 

i-.  vid:  c  .•  •  in  _  «!•.  »a~gi  =  dicut  miring  jc.nic  to  rcorre:ft  h»i  i^uiilsl 
••  it  :  .vii  iv  in  the  following  may: 


a'*  1  *•  '  .  .  s  a  cyclonic  :pt>i'ieu;tud  is  brought  up  l-y  (drip)  t  loud 

cc-d  d^rtcjned  j.n  the  coper  levels.  Iris  results  in  a 
,ruuoMvcj  ;•••.■>>  i  t  i  v  e  TCE  in  the  upper  layers  and  a  negative  TCE 
is:  lei  r  layers.  Such  processes  act  tJ  decrease  both  the  upper 
level  .  a  1 1  icy  r.  Ion; c  and  the  low  ]  t-<  el  cyclonic  -minds. 

•j  !  »  t.-liu-.v  of  T'-'E  (or  ;ur’  )  to  he  positive  uear  the  surface 

*  t> 

is  Ir.-.ely  due  to  downward  lotions  from  higher  layers  or  to 
*etf?iul  luji'uloat  mixing  processes  at  k-»cr  layers. 

C  ;  vritiiel  mixing  processes  aill  Cause  a  momentum  loss  at 

T:ic  level  of  luazircUja  cyclonic  wind  as  is  cib  jc  rvet 1  for  all  the 
li-  date.  sets. 

d )  .ti  tu-e-gii  up'gredient  uvacutUia  transport,  possibly  caused  by 
sr.sl  1  -scale  dynamic  effect,  might  occur  in  individual 

•- ■■  -  s,  steit  and  act  to  in-veiss  the  lovul  vertical  mir'd 
-i.  -  -  ;  us  ved  *»:.’d  suggested  Tor  individual  7  Sit. -use  Systems 

:  s  r  al_.  1V76;  Mullet  ruff  and  ‘biles',  1076),  the  r.ot 
ic  1  1  a.a’e  dynamic  effect  of  all  syssvcCtive  sys^rai  on  the 
t  'j  •  ^c  ,_ale  luc.acutuui  fields  is  uut  that  important  aud  the 

duv't-  gradient  mixing  process  appear*  to  be  dominant. 

5.11  Moan  „ni  hudy  Vertical  Transports 

lie  mean  (~  dtwllvl/op)  and  eddy  (—  elm'v']/6p  or  ouF  )  vertical 

W 

transports  of  Tangential  momentum  appear  to  exhibit  the  s&ibc  pattern  for 
all  t at  dev loping  cases,  with  the  stronger  systems  having  l«recr 
li'.gur.cuee  itvicr  to  Figs.  5—  6  and  Figs.  23  -.14).  To  compare  these  two 
transport  processes,  su  average  of  all  the  developing  systems  (WVD2, 
h’PD3 ,  ***!;■*  ( U j ,  A ’nil,  All'2.  Hlli3  and  ATL*4  (u)  )  has  been  taV.cn,  Results 
-re  shorn  in  Fig.  31.  ’Ihe  total  convective  effect  (or  the  sum  of  these 
two  p.-cce *•-«*?  is  also  shown.  Note  that  a  '  1  these  data  Sets  (which  are 


averaged  Together,' 


.deut.  This  figure  »bums  thc.t  Custtiui 
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Fig.  31.  Averaged  vertical  profiles  of  total  convective  effect  or  TCE 
(solid),  and  that  part  of  the  TCE  accomplished  by  the  Kean 
vertical  motion  or  -(3  [<*»]  tv] /dp)  (dashed)  and  that  part 
accomplished  by  the  eddy  process  or  .C"o«  (dotted)  for  all  the 
developing  data  sets  in  the  region  1-5  latitude  radius  around 
the  system  center.  The  heavy  dashed  curve  is  when  the  surface 
friction  is  added  to  the  cumulus  friction. 


vertical  transport  in  magnitude  and  its  pattern  is  quite  similar  to  the 
TCE.  These  results  indicate  that  cumulus  friction  is  an  important 
physical  process  and  ought  to  be  taken  into  account  in  numerical  model 
simulations  of  tropical  cyclones. 

The  pattern  of  this  cumulus  friction  agrees  well  to  that  of  the 
apparent  vorticity  source  calculated  by  Williams  and  Gray  (1973)  and 
Reed  and  Johnson  (1974).  Their  results  show  a  positive  momentum  or 
vorticity  source  at  the  upper  outflow  layer  and  a  negative  source  at 
middle  and  lower  layers  (surface  friction  was  not  taken  out  in  their 
calculation).  When  the  surface  friction  is  added  to  the  present  cumulus 


friction.  results  show  (heavy  dashed  curve  in  Fig.  31)  only  a  stronger 
low  level  negative  acceleration. 

Alternatively,  the  symmetric  equation  of  tangential  momentum  can  be 
*rit; in  et  advection  fora  as: 


Ihl 
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-here  =  ♦  “)  ♦  f  is  the  absolute  vorticity.  The  terra  —  “ 

'a  ■  pf  r  r  oe 

has  beer-  neglected  in  Eq.  26  because  it  is  always  much  smaller  than 
(.  .her  tc'r.s.  The  first  term  on  the  right  hand  side  represents  the 
horizc-utsl  advection;  and  the  second  and  the  third  terms  are  the 
vertical  advection  (VA)  and  vertical  eddy  transport  (or  cuF0) 
respectively.  Vertical  profiles  of  the  last  two  terms  are  shown  in  Fig. 
32.  Note  that  the  sum  of  these  two  terms  do  not  vertically  integrate  to 
zero.  Equation  26  provides  a  simpler  way  to  see  how  the  vertical  motion 
processes  can  act  to  produce  a  negative  tangential  acceleration  through 
a  deep  layer  -  up  to  300  to  350  mb.  A  similar  approach  has  been 
previously  taken  by  Gray  (1967)  to  study  momentum  changes  within  the 
eye-wall  regions  of  hurricane*. 

5.3  Environmental  Wind  Fields  and  Cumulus  Friction 

Because  cumulus  convection  is  primarily  a  down-gradient  mixing 
process,  the  total  convective  effect  mast  correspond  to  the 
i'.h i ror>  i;ntal  vertical  wind  shear  and  the  intensity  of  the  convection. 

3. ace  the  ne;>n  convective  effect  hat  almost  the  same  vertical  pattern 
{.-•  til  the  data  sets,  the  change  in  the  vertical  pattern  of  TCE  appeare 
-.o  he  largely  de  tcreined  by  the  change  in  the  cmsulus  friction  profile. 

’  ]  re,  if  the  vertical  mass  recycling  is  assumed  to  have  the  same 
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r.i  that  .shewn  in  Fig.  18  for  all  the  data  sets,  we  can  limit  our 
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Fig.  32.  Averaged  vertical  profiles  of  cumulus  friction  (dashed), 

vertical  advection  (VA)  tangential  momentum  by  mean  vertical  motion 
(dotted)  and  the  sum  of  these  two  (solid)  for  all 
the  developing  data  sets  in  the  region  1-5  latitude  radius 
arouud  the  system  center. 

discussion  to  the  relationship  between  the  cumulus  friction  profile  and 
the  mean  vertical  motion  and  tangential  wind  fields.  To  study  the 
dependence  of  the  cumulus  friction  on  the  environmental  wind  fields  data 
sets  are  grouped  into  3  different  intensity  stages.  Each  stage  includes 
both  developing  (D)  and  non-developing/steady  (N)  cases  as  shown  in 
Table  6.  Note  that  WPD4  and  ATD4  (or  III-N  category)  are  not  deepening 
but  have  a  steady  or  filling  tendency  trend. 
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The  data  sets  which  are 

TAHLE  6 

grouped  into  different 

Developing  (D) 

categories. 

Non-developing/steady  (N) 

(I)  Initial  Stage 

(I-D) 

(I-N) 

WPD2.ATD1 

WPN1.ATN1 

(II)  Intermediate  Stage 

(II-D) 

(II-N) 

WPD3 , ATD2 , ATD3 

WPN3.ATN3 

(III)  Mature  Stage 

(I II-D) 

(III-N) 

WPD4 (D) , ATD4(D) 

WPD4.  ATD4 

Figures  33  and  34  illustrate  the  averaged  vertical  profiles  (for  I, 
II  and  III  stages)  of  cumulus  friction,  tangential  wind,  and  mean 
veitical  motion  for  the  developing  (D)  and  non-developing/f il 1 ing  (N) 
cases.  The  profiles  of  each  of  these  three  parameters  share  similar 
vertical  patterns  for  the  D  and  N  cases,  though  the  magnitudes  are 
different.  These  cumulus  friction  profiles  are  similar  to  those 
diagnosed  hy  Hack  and  Schubert  (1980)  except  at  lower  levels  where  our 
developing  cases  have  a  smaller  amplitude. 

For  the  developing  cases,  the  cumulus  friction  does  not  change  its 
vertical  pattern  from  initial  to  mature  stage.  Its  amplitude  however 
increases  substantially  with  intensity  (Fig.  35a).  This  is  qnite 
reasonable  because  the  mean  vertical  motion  and  vertical  tangential  wind 
profiles  are  not  much  different  among  developing  systems  (Figs.  35b, 
35c),  and  thus  the  cumulus  friction  patterns  do  not  vary  much.  The 
increasing  magnitude  of  cumulus  friction  is  just  a  result  of  the 
increase  in  the  vertical  wind  shear  (see  Eq.  15). 

It  is  more  complicated  to  discuss  the  non-developing/steady  cases 
since  none  of  the  individual  composite  profiles  resembles  the  average 
profiles  of  Fig.  33.  Note  that  the  intensity  classifications  (I-N, 


Fig.  33.  Averaged  vertical  profiles  of  cumulus  friction  (cuF^)  for  all 
developing  (D)  and  non-developing/steady  (N)  cases  in  the 
region  1-5  latitude  radius  around  the  system  center. 
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Fig.  34.  Averaged  vertical  profiles  of  mean  vertical  motion  and  tangen¬ 
tial  wind  for  developing  (solid)  and  non-developing  (dashed) 
data  sets  in  the  region  1-5°  latitude  radius  around  the  system 
center. 
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Fig.  36.  1-3°  latitude  radial  and  azimuthal  averages  of  (a)  cumulus 

friction,  (b)  mean  vertical  motion  and  (c)  tangential  wind  for 
the  non-developing  tropical  cyclone  systems  of  cloud  cluster 
(I),  tropical  depression  or  storm  (II),  and  mature  typhoon  or 
hurricane  (III)  stages. 
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II- N  and  III-N)  are  not  a  progressive  increase  as  with  the  developing 
cases.  The  mean  vertical  motion  and  vertical  tangential  wind  profiles 
are  shown  in  Figs.  36b, c.  Obviously,  the  patterns  of  these  profiles 
differ  from  one  another  and  also  from  those  of  the  developing  cases. 

This  results  in  the  different  profiles  of  cumulus  friction.  For  the 
initial  stage,  not  much  difference  exists  in  the  mean  vertical  motion 
patterns  (except  the  magnitude)  between  the  I-N  and  I-D  cases;  however, 
the  tangential  wind  profiles  do  show  important  differences.  The  I-D 
case  (also  II-D  and  III-D)  has  the  level  of  maximum  tangential  wind  at 
about  850  mb  and  little  vertical  wind  shear  from  600-900  mb.  Contrary 
to  the  I-D  case,  the  tingential  wind  of  the  I-N  case  has  a  maximum  at 
about  500-700  mb  and  decreases  to  zero  (or  slightly  negative)  at  the 
surface.  The  positive  cumulus  friction  for  the  I-N  case  from  700  mb 
down  to  the  surface  is  probably  a  result  of  this  wind  shear  in  the 
middle  to  low  levels  and  the  effect  of  downdrafts  or  downward  mass  flux. 

For  the  intermediate  stage,  a  slightly  larger  600-900  mb  vertical 
wind  shear  is  found  for  the  II-N  compared  with  the  II-D  case.  The  most 
appreciable  difference  is  in  the  mean  vertical  motion.  Since  the 
vertical  mass  recycling  strongly  decreases  with  height,  the  decrease  of 
the  mean  vertical  motion  at  400-600  mb  for  the  II-N  case  implies  a 
secondary  divergent  layer  other  than  the  major  upper  outflow  layer. 

This  secondary  divergent  layer  and  the  vertical  wind  shear  yield  a 
positive  value  of  cumulus  friction  between  400-650  mb  for  the  II-N  case. 
For  the  III-N  case  the  patterns  of  tangential  wind,  vertical  motion,  and 
cumnlus  friction  are  quite  similar  to  those  of  the  developing  case 
(III-D),  But  the  magnitude  of  the  cumulus  friction  is  larger  for  the 

III- N  case.  This  is  a  result  of  the  greater  vertical  wind  shear  of  this 
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system. 

The  above  discussions  again  support  the  contention  that  cumulus 
convection  is  primarily  a  down-gradient  mixing  process; 'and  hence  the 
cumulus  friction  profile  has  to  correspond  to  the  environmental  flow 
features.  Therefore,  in  a  numerical  model,  cumulus  friction  ought  to  be 
determined  (or  parameterized)  according  to  the  simultaneously 
corresponding  flow  features  and  surface  energy  fluxes  (or  evaporation) 
manifested  through  the  vertical  mass  recycling  (as  discussed  in  Chapter 
4). 

5.4  The  Effects  of  Cumulus  Convection  on  the  Upper  Outflow  Layer 

The  net  effect  of  cumulus  convection  on  the  upper  outflow  layer  is 
to  decrease  the  anticyclonic  winds.  This  is  equivalent  to  an  increase 
in  absolute  angular  momentum  of  the  outward  moving  air  parcels.  An  air 
parcel  has  to  conserve  its  absolute  angular  momentum  along  its 
trajectory  unless  acted  upon  by  a  torque.  However,  trajectory  analyses 
for  hurricane  (ATD4)  and  typhoon  (WPD4)  systems  show  an  increase  of 
absolute  angular  momentum  associated  with  outflowing  air  parcels. 

The  trajectory  analysis  on  the  150  mb  isobaric  surface,  the  level 
of  maximum  concentrated  outflow,  has  been  carefully  made.  The  composite 
storm  is  assumed  to  be  in  a  steady  state  and  moving  with  constant 
direction  and  speed.  Streamline  patterns  are  obtained  at  two  different 
12  hour  times,  and  trajectories  constructed  (see  Saucier,  1955,  p.  313 
for  a  description  of  the  procedures  used).  Figure  37  depicts  an  example 
of  how  these  trajectories  are  determined.  Such  an  analysis  was 
performed  for  the  typhoon  (WPD4)  and  hurricane  (ATD4)  data  sets.  Figure 
38  shows  the  trajectories  and  the  constant  absolute  angular  momentum 
contours  at  150  mb  for  the  mean  typhoon  and  mean  hurricane.  Absolute 
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Storm  Center : 

P,  present  time 

P2  12  hours  later 
Streamlines  : 

dotted  curves  (present  time) 
dashed  curves  (12  hours  later) 


Fig.  37.  Schematic  diagram  showing  how  the  outflow  trajectory  was 
constructed.  Curves  ab  and  ad  are  the  paths  that  an  air 
parcel  will  travel  in  12  hours  from  point  a  if  it  follows  the 
stream  lines  of  the  present  time  (ab)  and  of  12  hours  later  (ad), 
respectively.  The  length  of  ab  and  ad  is  determined  by  the  wind 
speed.  The  curve  ac  is  the  trajectory  that  an  air  parcel  will 
travel  from  point  a  in  this  12  hour  period  with  respect  to 
the  moving  storm  center  P. 


angular  momentum,  N^,  is  defined  as: 

M  =  1/2  f  r2  +  rv 

a  o 

where  f  is  the  Coriolis  parameter  at  the  system  center.  These  figures 
o 

show  that  the  absolute  angular  momentum  of  an  air  parcel  generally 
increases  as  they  move  outward  from  the  storm  center. 

For  an  individual  air  parcel,  however,  there  may  still  be  a 
pressure  gradient  force  acting  upon  it.  This  will  also  cause  w  change 
of  absolute  angular  momentum  along  the  trajectory.  But  when  all  the 
trajectories  crossing  a  certain  radius  are  considered  the  net  tangential 
pressure  gradient  force  is  zero  and  the  absolute  angular  momentum  will 
still  increase  along  the  trajectory,  especially  at  inner  regions.  This 
suggests  the  importance  of  the  TCE  in  modifying  the  radial  circulation 
at  the  inner  region  (r  <  5°  latitude)  around  the  system  center.  It  is 
also  worthwhile  to  notice  that  the  greatest  crossing  angle  occurs  inside 
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about  5°  latitude  radius  and  to  the  south  of  the  storm  center. 

5.5  Discussion 

It  has  been  generally  accepted  that  the  low  level  convergence  in  a 
tropical  cyclone  system  is  primarily  due  to  surface  frictional  effects. 
However,  the  composite  data  show  that  for  a  developing  storm  the  inflow 
can  extend  up  to  300  mb.  This  deep  inflow  layer  is  probably  an 
important  factor  in  differentiating  a  developing  system  from  a  non¬ 
developing  one.  Figures  39  and  40  show  a  cross-section  of  the 
az imuthal ly-averaged  radial  wind  for  developing  vs.  non-developing 
systems  in  both  ocean  basins.  It  is  obvious  that  the  developing  systems 
have  a  much  deeper  inflow,  especially  within  the  inner  2-4°  radius 
region.  Note,  however,  that  the  non-developing  systems  have  a  more 
concentrated  low  level  inflow. 

It  is  possible  that  this  radial  circulation  (deep  inflow  layer)  is 
affected  by  the  cumulus  convection,  as  has  been  discussed  in  the  last 
section.  Eliassen  (1951)  showed  that  a  redial  circulation  has  to 
respond  to  an  angular  momentum  or  heat  source  in  a  stable  circular 
vortex.  The  pattern  of  this  radial  circulation  is  largely  determined  by 
the  inertial  stability  of  the  vortex  flow  (see  also  Schubert  and  Hack, 
1981).  Generally  speaking,  this  radial  circulation  can  be  considered  as 
a  secondary  circulation  to  wipe  out  the  imbalance  between  the  large 
scale  pressure  and  wind  fields.  The  present  results  again  show  similar 
patterns  between  the  radial  wind  and  the  TCE  profiles  (refer  to  Figs. 

27,  28,  39  and  40).  Therefore,  the  mutual  interaction  between  the  small 
scale  cumulus  convection  and  the  large  scale  radial  (and  tangential) 
flow  features  during  tropical  cyclone  development  appears  to  be  very 
important.  However,  the  exact  relationship  between  the  radial  wind  and 
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the  TCE  is  quite  complicated  and  this  observational  study  cannot  fully 
address  this  problem.  Special  numerical  modeling  work  along  these  lines 
would  be  highly  desirable. 

The  importance  of  this  radial  circulation  is  highlighted  when  the 
horizontal  transport  of  angular  momentum  is  considered.  As  has  been 
shown  in  Fig.  22a,  the  tangential  wind  increases  through  a  very  deep 
layer  at  the  inner  region  causing  the  increase  in  middle  to  upper  level 
vertical  wind  shear  as  the  system  intensifies.  This  increased  wind 
shear  will  then  be  able  to  sustain  a  greater  upper  level  warm  core,  and 
a  subsequent  lower  surface  pressure.  Since  cumulus  convection  is 
primarily  a  down  gradient  mixing  process  this  increase  of  vertical  wind 
shear  should  be  attributed  to  the  horizontal  transport.  The  importance 
of  the  horizontal  transport  of  angular  momentum  in  tropical  cyclone 
systems  has  also  been  discussed  by  Holland  (1981).  The  point  being 
addressed  here  is  the  importance  of  the  deep  inflow  and  the  concentrated 
outflow  layers  for  the  developing  cyclones.  This  deep  inflow  appears  to 
be  responsible  for  the  increase  of  tangential  wind  through  a  deep  layer. 
At  the  upper  levels,  the  outward  transport  of  tangential  momentum  just 
balances  the  large  increase  of  tangential  momentum  due  to  convective 
transport. 

Figure  41  illustrates  the  idealized  tangential  wind  changes  at  the 
inner  region  associated  with  the  horizontal  and  vertical  transports  (and 
surface  frictional  effect).  Curve  a  is  the  assumed  initial  tangential 
wind  profile  and  curve  b  is  the  resulting  tangential  wind  profile  after 
cumulus  convection  (or  local  mass  recycling)  if  we  assume  that  there  it 
no  net  radial  wind.  (The  differences  between  curves  a  and  b  at  each 
layer  have  been  scaled  to  represent  the  calculated  TCE.)  However,  there 
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Fig.  41.  Idealized  diagram  showing  the  effects  of  the  vertical  and  hor¬ 
izontal  transports  of  tangential  wind  and  of  the  surface  friction  on  the 
alteration  of  the  tangential  wind  profile  in  the  inner  region  of  a 
developing  tropical  cyclone  system.  The  bottom  diagram  is  the  combina¬ 
tion  of  all  the  processes  of  the  top  diagrams.  Curves  a  and  d  represent 
the  observed  initial  and  final  tangential  wind  profiles  for  a  tropical 
cyclone  system  which  is  undergoing  development  (top  right  diagram). 

Curve  b  is  the  speculated  profile  after  cumulus  convection  only  - 
(assume  no  radial  wind).  The  difference  between  curves  b  and  a  are 
scaled  to  represent  the  calculated  TCE  profile.  The  horizontal  tran¬ 
sport  of  tangential  wind  due  to  radial  circulation  will  change  the 
tangential  wind  profile  from  curve  b  to  curve  c.  At  the  same  time,  sur¬ 
face  frictional  effect  will  reduce  low  level  tangential  wind  from  curve 
c  to  curve  d.  Curve  d  is  thus  the  final  and  stronger  tangential  wind 
profile  resulting  from  all  these  effects. 
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is  always  a  mean  radial  circulation  accompanying  the  cumulus  convection. 
The  horizontal  transport  of  tangential  wind  associated  with  this  radial 
circulation  will  change  the  tangential  wind  profile  from  curve  b  to 
curve  c.  Along  with  the  vertical  and  horizontal  transport,  the  surface 
frictional  effect  always  acts  to  decrease  the  boundary  layer  tangential 
wind  from  curve  c  to  curve  d.  Curve  d  is  thus  the  final  and  stronger 
tangential  wind  profile  resulting  from  all  these  effects. 

At  the  outer  region,  the  increase  of  upper  level  anticyclonic  wind 
(Fig.  22b)  is  primarily  due  to  the  concentrated  outflow  (Figs.  39a  and 
40a)  which  exports  very  large  amount  of  absolute  vorticity  (which  is 
always  positive)  to  larger  radii  (Eq.  26).  In  the  middle  and  lower 
levels,  horizontal  and  vertical  transports  are  relatively  smaller  and 
tend  to  balance  each  other,  causing  only  a  slight  increase  in  tangential 
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6.  SUMMARY  AND  CONCLUSION 

Cumulus  friction  (or  the  unre solvable  part  of  the  convective 
effect)  can  be  calculated  as  a  residual  of  the  large  scale  momentum 
budget.  Such  calculations  have  been  performed  for  16  data  sets 
corresponding  to  different  intensity  stages  of  tropical  cyclones  in  the 
northwest  Atlantic  and  northwest  Pacific.  Azimuthal  averaging  has  been 
used  to  avoid  calculating  the  pressure  gradient  term.  To  reduce  hor¬ 
izontal  eddy  effects  on  the  residual  calculation,  only  the  1—5  latitude 
radii  results  have  been  analyzed.  These  residual  cumulus  friction 
effects  were  found  to  exhibit  the  same  pattern  for  all  the  developing 
data  sets:  positive  at  upper  and  lower  layers  and  negative  at  levels  in 
between.  The  stronger  systems  have  a  greater  amplitude  of  cumulus  fric¬ 
tion  than  the  weaker  systems.  The  differences  between  the  developing 
and  non-developing  systems  were  found  to  be  consistent  at  each  intensity 
stage  for  both  ocean  basins. 

Cumulus  friction  was  found  to  have  a  magnitude  comparable  to  the 
mean  vertical  motion  transport  and  a  pattern  similar  to  the  total  con¬ 
vective  effect.  Thus,  any  numerical  modeling  work  which  attempts  to 
simulate  tropical  convective  systems  ought  to  include  this  cumulus  fric¬ 
tion  effect.  To  properly  include  these  vertical  mixing  influences  in  a 
numerical  model  one  is  required  to  develop  an  adequate  parameterization 
scheme.  A  simple  single-cloud  model  was  developed  for  this  purpose. 

This  model  included  a  single  cloud  updraft  and  three  different  types  of 

i 

downward  motion.  The  results  show  that  a  simple  relationship  for  the 

influence  of  cumulus  convection,  namely  -  d[a>  (v  -v)]/dp,  can  be  used 

c  c 

to  represent  satisfactorily  such  cumulus  vertical  momentum  mixing  influ¬ 
ences.  In  addition,  the  consistency  between  the  calculated  residuals 
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and  model-diagnosed  results  add  support  to  both  observational  and 
model-inferred  results. 

The  total  convective  effect  (TCE)  was  found  to  exhibit  a  negative 
value  at  the  level  of  maximum  cyclonic  wind  and  a  positive  value  at  the 
primary  outflow  layer  (or  to  decrease  both  lower  level  cyclonic  and 
upper  level  anticyclonic  winds).  This  leads  to  the  conclusion  that 
although  up-gradient  transport  might  occur  under  certain  conditions  to 
increase  the  vertical  wind  shear,  the  net  influence  of  cumulus  convec¬ 
tion  in  tropical  cyclones  is  primarily  a  down-gradient  mixing  process 
which  acts  to  reduce  the  large  scale  vertical  wind  shear.  These  find¬ 
ings  are  consistent  with  the  recent  Harvard  School  speculations  on  the 
bulk  influences  of  cumulus  convection  on  large  scale  momentum  fields 
(Schneider  and  Lindzen,  1976;  Stevens  <5t  a_l. ,  1977;  Shapiro,  1978; 
Shapiro  and  Stevens,  1980). 

The  findings  of  this  paper  do  not  support  the  Imperial  College, 

O.K.  ideas  of  Moncrieff  and  Green  (1972)  and  Moncrieff  and  Miller  (1976) 
on  the  possible  influence  of  deep  convection  acting  to  increase  vertical 
wind  shear  (at  least  in  an  aggregate  sense)  or  to  the  hypothesis  of  Gray 
(1979)  and  McBride  (1979)  that  an  up-gradient  transport  of  momentum  by 
deep-cumulus  may  be  an  important  mechanism  in  tropical  cyclone  genesis. 
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LIST  OF  SYMBOLS 


Bracket  denotes  grid  box  average 

Prime  denotes  deviation  from  grid  box  average 

Overbar  denotes  azimntbal  mean 

Condensa  tion 

Drag  coefficient 

1.  Speed  of  cyclone  movement 

2.  As  a  subscript  denotes  cloud 

As  a  subscript  denotes  cloud  downdraft 
Cumulus  friction  (tangential  component) 

Fraction  of  total  downward  motion  in  form  of  nth  kind  of 
downdraft 

As  a  subscript  denotes  downward  motion 

Evaporation 

Surface  evaporation 

As  a  subscript  denotes  environmental  cloud  free  area 
Subgrid  scale  effect  (tangential  component) 

Coriolis  parameter 

Coriolis  parameter  at  the  storm  center 

1.  Gravity 

2.  As  a  subscript  denotes  environmental  gentle  subsidence 

As  a  subscript  denotes  vertical  levels,  k  =  0  at  the 
surface 

Absolute  angular  momentum 

As  a  subscript  denotes  mean  vertical  mass  flux 
As  a  subscript  denotes  meso-scale  downdrafts 
As  a  subscript  denotes  nth  kind  of  downdrafts 
Pressure,  vertical  axis 


Q 


Moisture  flux 
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LIST  OF  SYMBOLS  (cont’d) 

Mixing  ratio 
Satnration  mixing  ratio 

Condensation  to  rain  by  all  the  upward  motions 

1.  Radius,  radial  axis 

2.  As  a  subscript  denotes  local  mass  recycling 
Surface  friction  (tangential  component) 

Total  convective  effect  (tangential  component) 

1.  Radial  component  of  momentum,  radial  wind 

2.  As  a  subscript  denotes  upward  motion 

Tangential  component  of  momentum,  tangential  wind 
Sub-grid  scale  pressure  gradient  force 

1.  azimuth,  measured  counterclockwise  from  the  north 

2.  as  a  subscript  denotes  tangential  component 

Density 

Fraction  of  total  area  covered  by  cloud 
Stress  at  the  top  of  mixing  layer  (=0) 

Surface  stress 
Geopotential  height 

Vertical  p  velocity  or  vertical  mass  flux  (mb/d) 


Absolute  vorticity 
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